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EXECUTIVE  SUMMARY 

In  1989  a  study  was  initiated  to  document  the  geographic  distribution  and  relative 
abundance  of  crayfish  in  south-central  Ontario  lakes  and  to  establish  whether  or  not 
crayfish  can  be  used  to  evaluate  the  health  of  inland  lakes.  Several  crayfish  species 
that  are  found  in  these  lakes  are  known  to  be  intolerant  of  poor  water  quality.  This 
report  summarizes  the  results  of  a  survey  of  100  lakes  of  varying  size  and  water  quality. 
This  information  will  be  used  to  estimate  the  normal  range  of  crayfish  abundance  in 
relatively  unimpacted  lakes  in  south-central  Ontario.  Subsequently,  crayfish  catches 
from  presumably  impacted  lakes  will  be  compared  to  this  normal  range  in  order  to 
assess  the  biological  significance  of  human  activities. 

From  1989  through  1994,  crayfish  relative  abundances  were  estimated  using  baited 
minnow  traps.  Generally  54  traps  were  set  for  a  single  night  in  a  given  lake.  Several 
lakes  were  re-sampled  to  examine  the  repeatability  of  the  resultant  abundance 
estimates.  In  total,  7  crayfish  species  were  captured  (i.e.,  Cambarus  bartoni.  C. 
robustus.  Orconectes  immunis.  O.  obscurus.  O.  propinquus.  O.  rusticus  and  O.  virilis). 
Up  to  4  crayfish  species  were  caught  in  any  lake,  although  no  crayfish  were  captured  in 
29%  of  the  lakes.  By  contrast,  one  species  was  caught  in  23  lakes,  two  in  33  lakes,  3  in 
12  lakes,  and  4  species  in  3  lakes. 

We  sampled  representative  lakes  in  7  tertiary  watersheds,  although  the  majority  of  the 
lakes  were  from  4  watersheds.  Three  crayfish  species  were  found  in  all  4  of  these 


watersheds:  C.  bartoni.  O.  propinquus  and  O.  virilis.  In  addition,  C.  robustus  was 
collected  in  3  of  the  4  watersheds.  By  contrast,  O.  innmunis  and  O.  obscurus  were 
found  in  only  two  and  one  watersheds,  respectively.  Of  the  3  watersheds  that  were 
minimally  sampled  (i.e.,  4,  2  and  1  lakes,  respectively),  one  lake  also  yielded  O. 
rusticus.  an  introduced  species  like  O.  obscurus.  Because  the  biological  impacts  of 
these  two  introduced  species  are  largely  unknown,  more  lakes  from  these  watersheds 
should  be  sampled.  This  additional  information  will  provide  a  correction  factor  to  adjust 
the  normal  range  in  crayfish  abundance  for  situations  when  a  new  crayfish  species 
invades  a  lake. 

In  those  lakes  where  several  crayfish  species  co-occurred,  no  single  species  was 
always  the  most  abundant.  As  a  result,  any  crayfish  species,  including  an  introduced 
species,  may  be  the  most  abundant  species  in  a  lake.  This  variation  complicated  our 
efforts  to  predict  crayfish  abundance. 

In  those  lakes  where  crayfish  were  caught,  the  average  catch  per  night  (i.e.,  catch  per 
unit  effort  or  CUE)  ranged  from  0.01  to  12.8  crayfish  per  trap.  Generally,  male  crayfish 
were  caught  in  greater  numbers  than  females  for  all  species.  Reproductive,  or  Form  I 
males  were  usually  the  largest  crayfish  caught.  As  expected,  C.  robustus  was  the 
largest  species  caught,  whereas  O.  propinquus  was  the  smallest.  These  two  species 
are  the  largest  and  smallest  crayfish  found  in  Ontario. 

Surprisingly,  no  crayfish  were  caught  in  many  of  the  lakes  in  two  of  the  watersheds. 
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Forty-five  per  cent  of  the  lakes  that  were  surveyed  in  the  Cameron  Lake  drainage  (i.e., 
2HF)  supported  no  crayfish  populations.  Similarly,  42  per  cent  of  31  lakes  that  were 
surveyed  in  the  Moon  River  and  Go  Home  River  watershed  (i.e.,  2EB)  contained  no 
crayfish  populations.  These  results  differ  markedly  from  the  neighbouring  watersheds 
(i.e.,  the  Georgian  Bay  tributaries  [2EB]  to  the  north  and  the  Severn  River  watershed 
[2EC]  to  the  south)  where  only  4  and  20  per  cent  of  the  lakes  did  not  support  crayfish 
populations.  The  reasons  for  this  high  frequency  of  lakes  without  crayfish  are  unknown. 
However,  the  long-term  impacts  of  acid  rain  have  been  implicated  in  the  disappearance 
of  crayfish  from  another  lake  in  the  2HF  watershed. 

Crayfish  presence  -  absence  and  relative  abundance  (i.e.,  CUE)  were  examined  with 
respect  to  morphological  and  chemical  features  of  the  lakes.  Simple  graphs  and 
Pearson  and  Spearman  rank  correlations  suggested  curvilinear  relationships  between 
crayfish  relative  abundance  and  various  environmental  variables.  As  a  result,  we  also 
examined  the  correlations  from  multiple  regressions  using  squared  and  cubed  terms  to 
account  for  this  curvilinearity  (i.e.,  polynomial  regressions). 

Simple  correlations  between  total  crayfish  catch  and  the  environmental  variables 
indicated  no  strong  relationships.  However,  the  catch  of  each  species  appeared  to  be 
influenced  by  different  environmental  variables.  For  example,  C.  bartoni  and  O.  virilis 
were  only  abundant  in  the  high-elevation  lakes,  whereas  the  catches  of  O.  virilis  and  O. 
propinquus  were  positively  correlated  with  lake  maximum  depth.  The  polynomial 
regressions  indicated  that  conductivity  and  Mg  concentration  were  correlated  with  the 
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abundances  of  C.  bartoni.  Similarly  the  catch  of  O.  propinquus  was  correlated  with  Al 
and  NH4  concentrations,  whereas  O.  virilis  abundance  was  correlated  with  Na,  TKN  and 
DOC.  As  a  result,  a  separate  model  would  be  required  to  predict  the  abundance  of 
each  species. 

We  also  employed  canonical  correspondence  analysis  to  explore  multivariate  linkages 
between  crayfish  abundances  and  the  environmental  data.  However,  few  of  the 
environmental  parameters  were  correlated  with  the  multivariate  results.  The  largest 
correlations  resulted  from  polynomial  regressions  using  parameters  such  as  elevation, 
alkalinity,  DOC,  Mn,  NH4,  TKN,  pH  and  SO4.  In  combination  these  environmental 
variables  may  implicate  the  acid-rain  problem  that  historically  impacted  many  of  the 
lakes  in  south-central  Ontario. 

Crayfish  presence  -  absence  data  were  evaluated  by  comparing  the  average  values  of 
the  various  environmental  factors  for  those  lakes  where  a  species  was  present  with  the 
corresponding  averages  for  lakes  where  that  species  was  absent.  For  example,  C. 
bartoni  was  found  in  deep,  high  elevation  lakes  that  were  somewhat  acidic  with  low 
DOC  concentrations.  Orconectes  virlis  populations  were  found  in  lakes  with  higher  pH, 
higher  total  phosphorus  and  lower  Mn  concentrations.  By  contrast,  most  lakes  with  O. 
propinquus  were  at  lower  elevation  on  average,  with  higher  pH,  DOC,  total  phosphorus, 
and  shallower  maximum  depths  with  lower  concentrations  of  SO4. 

One  of  the  goals  of  this  study  was  to  characterize  the  normal  range  in  crayfish 
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abundances  with  respect  to  relatively  unimpacted  lakes  in  south-central  Ontario. 
Unfortunately,  this  study  failed  to  reveal  any  simple  patterns.  Different  combinations  of 
crayfish  species  were  found  in  the  different  watersheds,  including  two  introduced 
species  (O.  obscurus  and  O.  rusticus).  Moreover,  an  unusually  large  number  of  lakes 
in  several  watersheds  did  not  support  any  crayfish  species.  The  abundances  and 
occurrences  of  the  different  crayfish  species  were  correlated  with  different 
environmental  variables.  For  example,  lake  elevation  and  lake  maximum  depth  were 
important  parameters  for  several  species.  In  addition,  a  number  of  acid-rain  related 
variables  were  also  important.  One  explanation  for  these  findings  rests  in  the  common 
belief  that  many  of  these  lakes  have  recovered  chemically  from  historical  damage  from 
acid  rain.  Perhaps  crayfish  species  have  been  unable  to  recolonize  these  lakes  despite 
improvements  in  water  quality.  If  crayfish  are  to  be  used  as  biomonitors  of  the  health  of 
small,  inland  lakes  in  south-central  Ontario,  it  is  clear  that  more  information  and  further 
data  analysis  are  required  before  we  will  fully  understand  the  complex  patterns 
observed  in  this  study. 
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INTRODUCTION 

Bottom-dwelling  invertebrates  in  lakes  and  streams  are  routinely  monitored  because 
changes  in  the  composition  of  the  benthic  community  often  reflect  changes  in 
ecosystem  health  (Lauritsen  et  al.  1985,  Thornley  1985,  Barton  1989,  Reice  and 
Wohlenberg  1993).  In  contrast  with  other  types  of  biomonitors,  invertebrate  taxa  are 
useful  as  early-warning  indicators  because  of  their  small  size,  relatively  short  life  spans 
and  varying  sensitivity  to  changes  in  environmental  conditions  (Johnson  et  al.  1993, 
Rosenberg  and  Resh  1993). 

Comparative  studies  provide  a  basis  to  evaluate  changes  in  the  composition  of  the 
benthic  community  (Underwood  1991,  1993).  For  example,  long-term  studies  quantify 
year-to-year  (or  temporal)  variation  (Osenberg  et  al.  1994),  whereas  broad  geographic 
surveys  measure  spatial  variation  (e.g.,  Harvey  and  McArdle  1986,  Stephenson  et  al. 
1994).  When  the  results  of  temporal  and  spatial  studies  are  combined,  the  normal 
range  of  variation  can  be  estimated  (e.g.,  Thompson  1938,  Underwood  1991,  Yan  et  al. 
1996).  Once  the  normal  range  has  been  established,  we  can  build  models  that  predict 
biological  responses  to  stresses  arising  from  human  activities  (Schindler  1987, 
Underwood  1991). 

In  1988,  the  Ontario  Ministry  of  the  Environment  (MOE;  now  called  the  Ontario  Ministry 
of  Environment  and  Energy,  MOEE)  initiated  a  biomonitoring  project  to  assess  long- 
term  trends  in  benthic  invertebrate  communities  in  softwater  Precambrian  Shield  lakes 
(for  details  see  Reid  et  al.  1994a,  b).  This  project  represents  one  subcomponent  of  the 
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federal-provincial  Long  Range  Transport  of  Airborne  Pollutants  (LRTAP)  programme 
that  was  designed  to  identify  anticipated  changes  in  the  health  of  inland  aquatic 
ecosystems  associated  with  reductions  in  sulphur  emissions  (Shaw  et  al.  1992, 
Stephenson  et  al.  1994).  A  standard  set  of  survey  protocols  is  used  by  members  of  the 
LRTAP  programme  to  permit  comparisons  among  data  sets  collected  from  6  areas  in 
eastern  Canada  (e.g.,  see  Shaw  et  al.  1 995).  The  goals  of  this  programme  include  the 
following:  (1)  to  estimate  the  geographic  ranges  of  common  species  offish  and 
invertebrates;  (2)  to  compare  relative  abundances  of  fish  and  invertebrate  taxa  among 
lakes  and  among  regions;  and  (3),  to  detect  and  ultimately  predict  the  response  offish 
and  invertebrate  species  assemblages  to  anticipated  changes  in  water  chemistry. 

Many  species  of  benthic  invertebrates,  including  snails,  clams  and  aquatic  insects, 
perish  in  acidic  waters  (Dillon  et  al.  1984,  Okland  and  Okland  1986,  Mackie  1989, 
Stephenson  et  al.  1994).  Several  crayfish  species  are  also  extremely  sensitive  to  low 
pH  (Malley  1980,  Berrill  et  al.  1985,  France  1987,  1993,  Davies  1989).  Consequently, 
one  subcomponent  of  the  LRTAP  programme  was  designed  to  assess  trends  in 
crayfish  populations  over  time.  Fourteen  lakes  in  south-central  Ontario  were  surveyed 
each  year  to  examine  temporal  trends  in  crayfish  catches  (Reid  and  David  1990).  Over 
the  first  5  years  crayfish  catches  in  many  of  these  lakes  declined  (David  et  al.  1994), 
although  the  cause  of  this  decline  was  not  determined. 

To  complement  the  temporal  or  time-trend  data,  a  second  project  was  initiated  in  1989 
to  identify  the  spatial  distribution  of  crayfish  in  representative  lakes  in  south-central 
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Ontario.  This  spatial  survey  provides  a  snapshot  of  crayfish  distributions  that 
presumably  reflects  crayfish  dispersal  after  the  last  glacial  retreat,  as  modified  by  biotic 
(i.e.,  predation  and  competition;  Olsen  et  al.  1991 ,  DiDonato  and  Lodge  1993,  Hill  and 
Lodge  1994)  and  anthropogenic  factors  (i.e.,  man-induced;  Davies  1989,  Gunn  and 
Keller  1990,  France  1993).  Seven  native  species  of  crayfish  are  assumed  to  have 
colonized  Ontario  lakes  and  streams  from  glacial  refugia  in  central  and  eastern  areas  of 
the  United  States  (Crocker  and  Barr  1968,  Guiasu  et  al.  1996).  Two  additional  species 
have  been  introduced  into  Ontario  waterways;  presumably  by  fishermen  using  crayfish 
as  bait  (i.e.,  Orconectes  obscurus  and  O.  rusticus:  Crocker  and  Barr  1968,  Berrill  1978, 
Momot  1992). 

This  spatial  survey  is  based  on  mid-summer  catches  of  crayfish  in  baited  traps.  The 
resultant  data  are  tabulated  as:  (1)  the  occurrence,  catch  per  unit  effort  and  the 
proportion  of  males  and  females  that  were  caught;  and  (2),  the  size  distribution  of  the 
catch.  Physical  and  chemical  data  for  each  lake  are  used  to  model  or  explain  observed 
patterns  in  species  composition  and  relative  abundance  with  respect  to  lake 
morphology  and  water  chemistry.  Further  development  of  these  models  will  permit:  (1) 
the  prediction  of  crayfish  species  composition  and  relative  abundance  in  lakes  that 
have  not  been  surveyed;  (2)  the  assessment  of  impacts  associated  with  changes  in 
water  and  sediment  chemistry;  and  (3),  the  definition  of  unimpacted  or  "normal"  lakes  in 
south-central  Ontario. 


METHODS 
Study  Area 

Most  of  the  survey  lakes  are  located  in  south-central  Ontario  in  the  county  of  Haliburton 
or  in  the  districts  of  Nipissing,  Muskoka  and  Parry  Sound.  The  nearest  major  urban 
centres  are  Peterborough  (120  km)  and  Toronto  (200  km)  to  the  south,  and  North  Bay 
(130  km)  and  Sudbury  (200  km)  to  the  north  (Figure  1).  The  majority  of  the  lakes  are  in 
4  tertiary  watersheds  of  the  Great  Lakes  that  Cox  (1978)  has  identified  as:  the  Georgian 
Bay  tributaries  (2EA);  Moon  River  and  Go  Home  River  (2EB);  the  Severn  River  (2EC); 
and  the  Cameron  Lake  drainage  (2HF).  A  limited  number  of  lakes  (7)  were  sampled  in 
3  additional  watersheds;  the  Georgian  Bay  Island  Tributaries  (2CF,  the  Sudbury  area); 
the  Scugog  River  (2HG,  the  Peterborough  -  Lindsay  area);  and  the  Upper  Madawaska 
River  (2KD,  the  Algonquin  Park  area).  The  geographic  co-ordinates  for  each  lake  are 
tabulated  with  simple  descriptions  of  lake  morphology  and  lake  elevation  (Table  1). 


Climate  and  Vegetation 

The  average  annual  precipitation  in  this  part  of  south-central  Ontario  is  90-1 10  cm  (i.e. 
a  30-year  average).  A  total  of  240-300  cm  of  snow  falls  each  year,  generally  between 
December  1  and  April  10.  January  temperatures  average  -10°  C,  whereas  the  July 
average  is  17.5°  C.  Mean  annual  temperature  is  about  5°  C.  The  lakes  are  generally 
frozen  from  the  first  week  of  December  to  mid-April  (Shaw  et  al.  1992). 


The  study  area  is  situated  in  the  Great  Lakes  -  St.  Lawrence  forest,  a  region 
characterized  by  white  pine  (Pinus  stEQbus),  red  pine  (P.  resinosa),  eastern  hemlock 
(Tsuga  canadensis)  and  yellow  birch  (Betula  alleghaniensis).  Other  common  species 
include  beech  (Fagus  grandifolia).  white  oak  (Quercus  alba),  sugar  maple  (Acer 
saccharum).  basswood  (Tilia  americana).  eastern  white  cedar  (T.  occidentalis),  red 
maple  (A.  rubrum).  red  oak  (Q.  rubra),  white  birch  (B.  papyrifera)  and  trembling  aspen 
(Populus  tremuloides).  Generally,  forest  stands  on  the  north-facing  slopes  are 
coniferous,  whereas  south-facing  slopes  are  characterized  by  deciduous  trees. 
Coniferous  trees  are  also  common  along  the  shore  of  most  lakes  (Reid  et  at.  1994a). 

Geology 

Most  of  the  study  lakes  are  situated  on  the  Canadian  Shield  (Jeffries  and  Snyder  1983, 
Girard  and  Reid  1990).  The  last  period  of  mountain  building  of  the  Precambrian  period 
was  the  Grenville  orogeny  which  formed  the  Grenville  province,  a  roughly  rectangular 
strip  of  land  400  km  wide  and  2000  km  long  stretching  from  Georgian  Bay  to  Labrador. 
The  bedrock  of  the  Grenville  province  is  granitic  in  composition,  usually  composed  of 
granitic  gneisses  and  migmatites  with  marble,  quartzite,  amphibolite  and  various 
igneous  intrusives  (pegmatites,  diorites,  and  metabasics)  of  less  importance  (Jeffries 
and  Snyder  1983).  During  the  Pleistocene  period  four  glacial  advances,  including  the 
most  recent  or  Wisconsin  glaciation,  have  rounded  and  polished  the  rock  outcrops  of 
the  Grenville  province,  lowering  the  relief,  scouring  many  river  and  lake  basins  and 
redepositing  glacial  debris,  including  moraines,  eskers,  outwash  sands,  gravel  and  lake 
sediments.  The  surficial  deposits  and  characteristics  of  the  watersheds  are  a 
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consequence  of  this  glacial  activity.  Minor  till  plains  and  thin  till  deposits  interrupted  by 
rock  ridges  are  the  dominant  surficial  features  (Jeffries  and  Snyder  1983,  Girard  et  al. 
1985,  Reid  and  Snyder  1986,  Girard  and  Reid  1990). 

Lake  Selection 

Four  tertiary  watersheds  were  selected  because  they  contain  lakes  that  are  sampled  in 
the  MOEE  time-trend  crayfish  study  (Reid  and  David  1990,  David  et  al.  1994).  Three 
other  watersheds  were  minimally  sampled:  2CF  (Lohi,  Clearwater,  Middle  and  Hannah 
lakes  near  Sudbury);  2KD  (Timberwolf  and  Louisa  lakes  in  Algonquin  Park);  and  2HG 
(Scugog  Lake  near  Peterborough).  Initially  candidate  lakes  were  selected  to  represent 
a  minimum  of  5%  of  the  total  number  of  lakes  within  a  size  range  of  10  to  99  hectares 
for  a  given  watershed  (using  counts  in  Cox  1978).  Final  selection  was  based  on 
accessibility,  which  was  generally  by  road.  This  quota  was  satisfied  for  most 
watersheds  by  the  end  of  the  1992.  Although  no  lakes  were  sampled  in  1993,  a 
number  of  larger  lakes  were  sampled  in  1994  (see  Table  1). 

Site  Selection 

After  arriving  at  a  lake,  the  field  crew  circumnavigated  the  lake  and  identified 
representative  littoral-zone  areas  that  characterized  the  dominant  substrate  types  (for 
details  see  Reid  et  al.  1994a).  Typically  the  nearshore  substrates  included  bedrock 
(BDR),  boulders  (BOU),  cobble  (COB),  pebble  (PEB),  silt  (SLT),  coarse  sand  (CSD) 
and  detritus  (i.e.,  organic  material  including  woody  debris:  LOS).  Rocks,  woody  debris 


and  macrophytes  provide  crayfish  with  refuges  from  predatory  fish  (Stein  1977,  Lodge 
et  al.  1986).  When  compared  with  other  substrate  types,  these  areas  tend  to  support 
greater  numbers  of  crayfish  whenever  predatory  fish  species  are  also  found  in  the  lake 
(e.g.,  Collins  et  al.  1983;  Kershnerand  Lodge  1995). 

After  circumnavigating  the  lake,  3  sites  were  selected  as  representative  areas  for 
subsequent  trapping.  Generally,  the  3  sites  were  characterized  as:  (1)  rocks  (and/or 
cobbles);  (2)  macrophytes;  and  (3),  detritus  including  logs.  Nearshore  substrates  in 
these  lakes  were  not  surveyed  by  SCUBA  divers  as  in  Reid  and  David  (1990).  As  a 
result,  detailed  descriptions  of  the  substrates  at  each  site  are  not  available. 

Water  Chemistry  Sampling  Procedures 

In  most  lakes,  water  samples  were  collected  on  the  same  day  that  the  traps  were  set  for 
crayfish.  However,  the  water  chemistry  of  several  of  the  lakes  is  regularly  monitored 
throughout  the  ice-free  season.  For  these  latter  lakes,  we  report  the  water  chemistry  as 
mean  values  for  the  ice-free  season.  These  ice-free  means  are  more  representative 
than  data  based  on  a  single  mid-summer  sample. 

Mid-lake  water  samples  were  collected  as  0-5  m  tube  composites  (see  Girard  and  Reid 
1990);  alternatively  3  L  of  water  were  collected  with  a  peristaltic  pump  (Masterflex 
Portable  Sampling  Pump,  Model  7570)  using  1/2  inch  i.d.  PVC  tubing  (Tygon  R3603). 
Depth  composites  were  obtained  by  pumping  water  into  an  8-L  polyethylene  carboy 
while  the  end  of  the  tube  was  lowered  to  the  desired  depth  and  brought  back  to  the 
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surface.  When  a  pump  was  used,  the  tube  was  lowered  and  raised  through  the  water 
column  at  a  constant  rate  to  sample  each  depth  equally.  The  resultant  water  was 
filtered  through  a  80-um  mesh  screen  placed  at  the  mouth  of  the  carboy  in  order  to 
remove  fine  particulates.  Subsamples  from  the  carboy  were  then  poured  into  the 
appropriate  bottles  and  containers  for  subsequent  analysis  at  the  laboratory.  The 
chemical  parameters  included:  alkalinity  (Alkti),  aluminum,  calcium,  chloride, 
conductivity  (Cond25),  dissolved  inorganic  carbon  (DIG),  dissolved  organic  carbon 
(DOC),  iron,  potassium,  magnesium,  sodium,  nitrogen  (NH4,  NO3  and  total  Kjeldahl 
nitrogen  [TKN]),  pH,  sulphate  (SOJ  and  total  phosphorus  (TP;  see  Table  2).  All 
chemical  analyses  were  completed  within  24  hours  of  sample  collection,  using  standard 
analytical  procedures  (see  Girard  and  Reid  1990). 

Fish  Community  Data 

The  Ontario  Ministry  of  Natural  Resources  (MNR)  in  Bracebridge,  Minden,  Parry  Sound 
and  Lindsay  provided  data  summarizing  the  fish  species  collected  in  surveys  of  many  of 
the  study  lakes  (Tables  3  and  4).  Fish  were  collected  using  various  types  of  gear 
including  experimental  gill  nets,  trapnets,  Windermere  traps,  minnow  traps  and  seine 
nets.  Because  several  common  fish  species  (principally  centrarchids  and  percids)  prey 
on  crayfish,  fish  inventories  are  helpful  when  interpreting  crayfish  catch  data.  Predatory 
fish  species  are  believed  to  reduce  baited-trap  CUE  estimates  because  crayfish  are 
less  active  in  lakes  where  predatory  fish  are  common  (Stein  1977,  Collins  et  al.  1983, 
Somers  and  Green  1993).  To  date,  no  correction  factors  are  available  to  adjust 


crayfish  CUE  data  to  compensate  for  the  presence  of  various  predatory  fish  species. 
Speculation  based  on  patterns  in  the  time-trend  crayfish  data  (David  et  al.  1994) 
suggests  that  year-to-year  variation  in  the  abundances  of  the  larger  predatory  fish  also 
affects  crayfish  CUE.  As  a  result,  the  future  development  of  any  correction  factors 
should  include  estimates  of  the  year-class  abundance  of  the  various  predatory  fish 
species. 

Crayfish  Sampling  Procedures 

We  used  the  same  trapping  methods  that  were  employed  in  the  MOEE  time-trend  lakes 
(Reid  and  David  1990,  David  et  al.  1994).  Crayfish  were  collected  with  standard, 
commercially  available,  wire-mesh  minnow  traps.  Funnel  entrances  were  enlarged  to 
3.5  cm  to  accommodate  large  crayfish  (Collins  et  al.  1983).  Traps  were  set  in  the 
afternoon  and  retrieved  the  following  day.  To  reduce  possible  bias  arising  from  crayfish 
escaping  during  the  trapping  period,  traps  were  set  and  picked  up  from  the  sites  in 
random  order. 

For  most  lakes  a  total  of  54  traps  was  set  for  a  single  night,  although  occasionally  some 
traps  were  not  found  the  next  day  (Table  5).  As  noted  above,  3  sites  were  selected  in 
each  lake  and  18  traps  were  set  at  each  of  the  3  sites.  From  1989  to  1992,  3  traplines 
were  set  perpendicular  to  the  shore  in  each  of  the  3  sites.  Successive  traplines  were 
set  approximately  5  m  apart.  Each  trapline  consisted  of  6  traps  set  at  1-m  depth 
intervals  from  1  to  6  m.  Traps  were  individually  marked  with  2-litre  plastic  bottles 
attached  to  a  measured  length  of  rope.  In  bays  that  were  less  than  6-m  deep,  the 
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remaining  traps  were  set  at  shallower  depths  to  ensure  a  constant  total  effort. 

Commencing  in  1994,  6  traps  were  individually  attached  to  a  nylon  line  at  3-m  intervals. 
Three  of  these  traplines  were  then  set  perpendicular  to  the  shore  for  a  total  of  18  traps 
per  site,  and  54  traps  per  lake.  This  design  is  similar  to  that  of  previous  years,  with  the 
exception  that  individual  depths  were  not  measured.  As  a  result,  trap  depths  ranged 
from  less  than  1  to  8  m  depending  on  the  slope  of  the  littoral  zone.  In  some  lakes,  the 
deepest  traps  were  in  shallower  water  than  in  earlier  years.  Consequently,  species 
such  as  Cambarus  bartoni  that  is  often  found  in  deeper  water  may  be  under- 
represented  in  the  1994  catches. 

Each  trap  was  baited  with  a  single,  perforated  plastic  film  canister  that  was  filled  with 
fish-flavoured  canned  cat  food  (David  et  al.  1994).  The  same  type  of  bait  was  used 
throughout  the  study  because  different  types  of  bait  will  affect  the  size  composition  and 
the  crayfish  species  that  are  caught  in  baited  traps  (Somers  and  Stechey  1986). 

On  each  sampling  period,  we  recorded  lake  name,  date,  site  and  trap  number  or  depth. 
The  species,  sex  and  carapace  length  of  each  crayfish  were  also  recorded.  Crayfish 
were  identified  using  taxonomic  keys  in  Crocker  and  Barr  (1968).  In  the  last  few  years, 
representative  specimens  were  donated  to  the  collections  at  the  Royal  Ontario 
Museum.  The  identification  of  these  crayfish  was  confirmed  by  Dr.  D.W.  Barr.  The  sex 
of  each  crayfish  was  classified  as  male  Form  I  (sexually  functional)  or  Form  II  (sexually 
non-functional),  or  female.  Dorsal  carapace  length  was  measured  with  vernier  calipers 
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in  mm  to  one  decimal  place.  Carapace  length  was  defined  as  the  distance  between  the 
tip  of  the  acumen  (the  extreme  anterior  end  of  the  rostrum)  and  the  posterior  end  of  the 
cephalothorax  (see  Illustration  in  Reid  and  David  1990). 

Baited  traps  were  set  in  approximately  20  lakes  each  year  from  1989  through  1994, 
although  no  "spatial"  lakes  were  surveyed  in  1993  (Table  5).  David  et  al.  (1994)  report 
trap-catch  data  for  14  time-trend  lakes  that  were  repeatedly  sampled  over  the  same 
years.  Because  the  catches  of  several  crayfish  species  declined  in  some  of  the  lakes  in 
that  time-trend  study,  7  of  the  100  lakes  were  sampled  twice  to  provide  a  modest  quality 
assurance/quality  control  assessment  (i.e.,  QA/QC).  Two  lakes  that  were  sampled  in 
1989,  1  in  1990,  and  4  lakes  that  were  sampled  in  1991  were  re-sampled  in  1994 
(Table  5).  These  lakes  were:  Bigwind,  Chub,  Dickie,  Moot,  Red  Chalk  East,  Red  Chalk 
Main  and  Twelve  Mile.  Catch  data  for  the  second  year  were  plotted  against  the  CUE 
from  the  first  year  and  the  results  were  compared  to  the  1 :1  line  to  evaluate  potential 
catch  biases  over  time. 

Data  Analysis 

Trap  catches  for  each  year  were  tabulated  by  species  and  by  sex  to  search  for  temporal 
trends  attributable  to  trapping  biases.  Similarly,  minimum,  maximum  and  mean 
carapace  lengths  of  the  catches  of  each  species  were  tabulated  by  reproductive  group 
for  each  lake.  Unfortunately,  there  were  often  very  few  animals  in  most  categories.  As 
a  result,  the  mean  values  and  ranges  are  probably  unrepresentative  for  most  lakes 
because  of  the  small  sample  sizes.  Where  more  than  10  animals  were  available  in  a 
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given  category,  we  correlated  mean  carapace  length  with  the  corresponding  water 
chemistry  variables  to  search  for  relationships  between  average  crayfish  size  and  water 
quality  (e.g.,  France  1985). 

The  CUE  of  the  most  frequently  encountered  species  (i.e.,  C.  bartoni.  C.  robustus.  O. 
propinquus  and  O.  virilis)  and  total  crayfish  CUE  were  also  correlated  with  lake 
morphology  and  water  chemistry  parameters  to  search  for  correlations  that  might 
indicate  potential  cause-and-effect  relationships.  Lakes  where  a  species  was  not 
encountered  were  excluded  from  this  analysis  (i.e.,  lakes  with  a  CUE  of  0  were 
excluded).  Both  Pearson  and  Spearman  rank  correlations  were  calculated  because 
linear  or  monotonic  (i.e.,  rank-order)  relationships  may  exist.  Multiple  regressions  were 
also  calculated  using  the  observed,  squared  and  cubed  forms  of  the  environmental 
variables  to  search  for  nonlinear  trends  in  the  CUE  (i.e.,  polynomial  regressions; 
MINITAB;  Ryan  et  al.  1988).  Traditional  multiple  regressions  using  all  of  the 
environmental  variables  simulataneously  were  not  utilized  because  of  problems  (e.g., 
instability)  due  to  multicollinearity.  To  facilitate  comparisons  between  the  Pearson, 
Spearman  and  polynomial-regression  results,  the  square-root  of  the  multiple  correlation 
was  tabulated  from  the  regression  analysis.  Similarities  among  the  3  correlation 
coefficients  were  used  to  identify  potentially  important  environmental  variables,  as  well 
as  the  nature  of  the  underlying  relationship  (i.e.,  linear,  curvilinear  or  nonlinear). 
Traditional  statistical  criteria  (e.g.,  P<0.05)  were  not  utilized  because  the  relatively  large 
number  of  correlations  are  expected  to  suggest  many  statistically  significant 
relationships  purely  by  chance.  Because  this  is  an  exploratory  analysis,  we  focus  on 
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correlations  that  exceeded  the  arbitrary  threshold  of  0.30  rather  than  statistically 
significant  correlations.  Those  environmental  variables  with  correlations  exceeding  this 
threshold  minimally  account  for  (or  explain)  approximately  10%  of  the  variation  in  the 
associated  crayfish  CUE. 

To  further  explore  the  crayfish-environment  relationships,  a  type  of  canonical 
correspondence  analysis  was  utilized  (e.g.,  CANOCO;  ter  Braak  1988,  1990,  Palmer 
1993).  To  this  end,  standard  correspondence  analysis  (i.e.,  a  CA  ordination)  was  used 
to  summarize  the  multivariate  trends  in  crayfish  CUE  and  the  associated  presence- 
absence  matrices  (NTSYS;  Rohlf  1993).  Correspondence  analysis  uses  a  chi-squared 
measure  of  inter-lake  and  inter-species  distance  which  is  consistent  with  the  expected 
curvilinear  distribution  of  species  abundances  along  environmental  gradients  (i.e., 
Gaussian  or  normal  curves;  ter  Braak  1988,  ter  Braak  and  Prentice  1988).  Because  of 
several  unusually  high  crayfish  catches,  CUE  values  were  down-weighted  by 
logarithmically  transforming  the  data.  Similarly,  the  undue  influence  of  rare  species  was 
down-weighted  by  adding  a  constant  (i.e.,  1.0)  to  each  value  in  the  log-transformed 
CUE  and  presence-absence  matrices  (cf.  ter  Braak  and  Prentice  1988,  Reid  et  al. 
1995).  Biplots  of  the  ordination  scores  for  the  6  crayfish  species  caught  in  more  than 
one  lake  (i.e.,  O.  rusticus  was  excluded)  and  the  71  lakes  where  crayfish  were  captured 
were  used  to  interpret  the  multivariate  patterns. 

The  first  3  CA  axes  were  correlated  with  the  combined  matrix  of  22  morphology  and 
water-chemistry  variables  using  Pearson,  Spearman  and  polynomial  regression 
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correlations  (MINITAB;  Ryan  et  al.  1988).  The  use  of  multiple  regression  with  CA  is 
similar  to  traditional  canonical  correspondence  analysis  (e.g.,  see  ter  Braak  1988, 
Palmer  1993).  Because  the  regression  analysis  eliminates  lakes  with  missing  data,  a 
small  number  of  missing  values  in  the  matrix  of  lake  morphology  and  water  chemistry 
variables  (e.g.,  see  Tables  1  and  2)  were  replaced  by  estimates  based  on  multiple 
regressions  using  all  of  the  remaining  environmental  variables  as  predictors. 

The  Pearson,  Spearman  and  polynomial  regression  models  were  used  to  examine 
relationships  between  crayfish  CUE  (and  the  multivariate  ordination  summaries)  and 
the  environmental  variables.  To  complement  the  correlation  and  regression  analyses, 
we  also  used  one-way  analysis  of  variance  (ANOVA)  to  compare  the  means  of  the 
environmental  variables  for  lakes  with  and  lakes  without  the  different  crayfish  species 
(MINITAB;  Ryan  et  al.  1988).  This  analysis  focuses  on  the  presence-absence  aspect  of 
the  crayfish  catches  in  an  attempt  to  establish  which  environmental  variables  best 
predict  whether  or  not  crayfish  will  be  found  in  a  given  lake. 


RESULTS  AND  DISCUSSION 

This  spatial  survey  provides  no  historical  information  as  to  whether  or  not  crayfish 
species  assemblages  have  changed  (cf..  Lodge  et  al.  1986,  France  and  Collins  1993). 
However,  Reid  and  David  (1990)  and  David  et  al.  (1994)  document  temporal  changes  in 
crayfish  abundances  in  a  related  study  of  14  south-central  Ontario  lakes  that  may 
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reflect  acid-rain  or  related  impacts  (i.e.,  declines  in  crayfish  abundance  were  observed 
in  several  lakes).  The  magnitude  of  the  year-to-year  changes  reported  by  David  et  al. 
(1994)  should  be  considered  when  evaluating  the  spatial  patterns  described  in  this 
study  because  this  survey  spanned  several  years. 

Elsewhere,  year-to-year  changes  in  crayfish  abundances  have  been  observed  in 
response  to  water-quality  changes  (e.g.,  France  1987,  Davies  1989,  Gunn  and  Keller 
1990)  and  the  invasion  of  introduced  species  (e.g.,  Capelii  1982,  DiDonato  and  Lodge 
1993).  Moreover,  France  and  Collins  (1993)  have  attributed  the  disappearance  of  O. 
virilis  from  Plastic  Lake  (one  of  the  temporal  reference  lakes  in  David  et  al.  [1994])  to 
changes  associated  with  a  gradual  decline  in  lake  pH  (see  Dillon  et  al.  1987).  Plastic 
Lake  has  not  been  recolonized  by  crayfish  since  the  last  crayfish  was  captured  in  1984 
(France  and  Collins  1993,  David  et  al.  1994).  The  factors  affecting  the  recolonization  of 
small  lakes  by  crayfish  have  not  been  identified;  however,  the  formation  of  ice  on  the 
bottom  of  the  inter-connecting  streams  during  the  winter  may  be  an  important  barrier  in 
this  part  of  Ontario.  For  example,  Aiken  (1968)  reported  that  O.  virilis  died  when  the 
substrates  of  small  Alberta  streams  froze. 

Lake  Morphology  and  Water  Chemistry 

Crayfish  catches  are  tabulated  by  tertiary  watershed  and  reported  with  information  on 
lake  morphology  (Table  1),  water  chemistry  (Table  2)  and  the  fish  communities  (Tables 
3  and  4).  Because  many  of  these  lakes  may  have  been  historically  impacted  by  acid 
precipitation  (Dillon  et  al.  1984,  Gunn  et  al.  1988),  acid  sensitivity  ratings  are  included 
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(Table  1).  The  majority  of  the  study  lakes  fall  within  the  moderate  sensitivity  (40  to  199 
A^eq-L'^)  or  the  extremely  sensitive  categories  (0  to  39.9  ^eq-f^).  As  a  result  of 
emission  controls,  the  water  chemistry  of  some  lakes  (especially  those  lakes  classified 
as  acidic  or  extremely  sensitive)  may  have  changed  over  the  last  several  decades 
(Dillon  et  al.  1987,  Gunn  and  Keller  1990).  However,  crayfish  populations  may  not  have 
recovered  in  (or  possibly  recolonized)  lakes  where  water  quality  has  improved  (e.g., 
France  and  Collins  1993). 

Watershed  2CF  -  Georgian  Bay  Island  Tributaries 

Lakes  in  watershed  2CF  drain  into  Georgian  Bay  and  Lake  Huron.  These  lakes  are 
unique  in  relation  to  many  of  the  other  survey  lakes  because  of  their  proximity  to  the 
Sudbury  smelting  industries  (e.g.,  see  the  Sudbury  Environmental  Study  1982).  Three 
of  the  4  lakes  that  were  surveyed  (i.e.,  Lohi,  Hannah  and  Middle)  were  treated  in  whole- 
lake  manipulations  (i.e.,  by  neutralization  and  fertilization)  in  the  1970s  (Yan  et  al. 
1996).  Prior  to  the  manipulations,  all  4  lakes  were  acidic  (i.e.,  pH  <  4.6;  Sudbury 
Environmental  Study  1982).  The  biological  recovery  of  these  lakes  is  ongoing  (Gunn 
and  Keller  1990,  Yan  et  al.  1996). 

The  mean  area  of  the  4  lakes  was  43.1  ha,  ranging  from  28.2  ha  to  76.5  ha  (Table  1). 
Mean  depth  ranged  from  4.0  to  8.4  m.  In  contrast  to  the  other  watersheds,  there  was 
little  difference  in  the  elevation  of  these  lakes  (range  from  303  m  above  sea  level 
[mASL]  to  307  mASL).  The  pH  values  ranged  from  5.37  to  7.23,  and  alkalinity  (i.e., 
total  inflection-point  alkalinity)  ranged  from  1.48  to  14.88  mg-L""  (Table  2).  DOC 
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concentrations  ranged  from  2.0  to  3.5  mg-L"^ ;  TKN  values  ranged  from  180  to  300 
A^g-L  ^  and  total  phosphorus  ranged  from  5.37  to  7.23  /^g-L"V  Aluminum  concentrations 
ranged  from  0  to  41  /ug-L\  whereas  sulphate  values  ranged  from  15.1  to  32.0  mg-L"\ 

Two  of  the  lakes  (i.e.,  Hannah  and  Middle)  had  the  highest  chloride  concentrations  of 
any  lakes  in  this  survey  (77.9  and  54.7  mg-L^).  These  two  lakes  also  had  the  highest 
sodium  values  (44.8  and  29.2  mg-L"^)  with  the  exception  of  Richmond  Lake  in 
watershed  2EA  with  a  sodium  concentration  of  50.2  mgL\  All  of  these  lakes  are  near 
major  highways.  As  a  result,  the  high  sodium  and  chloride  concentrations  probably 
reflect  runoff  of  road  salt  used  in  the  winter.  For  example,  95  to  98%  of  the  total  input 
of  Cr  to  Middle  and  Hannah  lakes  in  the  1970s  is  believed  to  have  originated  from  road 
salt  (Sudbury  Environmental  Study  1982). 

Watershed  2EA  -  Georgian  Bay  Tributaries 

Lakes  in  this  watershed  also  drain  into  Georgian  Bay  and  Lake  Huron.  Here  we 
sampled  5.4%  of  the  lakes  that  ranged  from  10  to  99  hectares  in  size  (i.e.,  22  lakes  of  a 
total  of  407;  Cox  1978).  The  mean  area  of  these  lakes  was  48.8  hectares  (ranging  from 
2.8  to  96.5  ha;  Table  1 ).  Because  some  of  the  lakes  in  this  size  range  were 
inaccessible  by  road,  5  larger  lakes  (138.9  to  363.4  ha)  were  also  sampled  (i.e.,  Harris, 
Rock  Island,  Fowke,  Manitouwaba  and  Grass  lakes).  The  mean  depth  of  all  27  lakes 
ranged  from  1.7  to  14.4  m.  Elevation  ranged  from  175  to  389  mASL. 

The  pH  of  these  lakes  ranged  from  4.92  (Black  Lake)  to  7.56  (Fowke  Lake),  whereas 
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alkalinity  ranged  from  -0.62  (Black)  to  18.38  mg-L'  (Grundy  Lake;  see  Table  2). 
Concentrations  of  DOC  ranged  from  3.0  (in  Black  and  Himbury  lakes)  to  8.3  mg-L"^  (in 
Grundy  and  Crawford  lakes)  and  TKN  ranged  from  180  (Kernick)  to  450  Aig-L'^  (Grundy 
Lake).  Total  phosphorus  ranged  from  3.0  (Manitouwaba  Lake)  to  13.85  ^g-L"^ 
(Whalley  Lake)  and  aluminum  concentrations  ranged  from  0  (Bells  and  Limestone 
lakes)  to  105  /^g-L"^  in  Round  Lake.  Sulphate  (SO4)  ranged  from  4.85  (Long  Lake)  to 
7.75  mg-L'^  (Fowke  Lake). 

Watershed  2EB  -  Moon  River  and  Go  Home  River 

Lakes  in  this  watershed  drain  into  Georgian  Bay  and  Lake  Huron  through  the  Moon  and 
Go  Home  rivers.  We  sampled  3.8%  of  the  lakes  ranging  from  10  to  99  ha  in  this 
watershed  (i.e.,  21  of  a  total  of  556;  Cox  1978).  The  average  area  of  these  lakes  was 
43.8  ha  with  a  range  of  8.7  to  93.6  ha  (Table  1).  Less  than  5%  of  the  lakes  in  the  10-to- 
99  ha  range  were  sampled  because  of  limited  road  access.  As  a  result,  1 0  larger  lakes 
(111  to  763.4  ha)  were  also  sampled.  The  average  mean  depth  of  all  31  lakes  was  5.0 
m,  with  a  range  of  2.2  to  10.7  m.  Lake  elevation  ranged  from  194  to  401  mASL. 

The  pH  of  these  lakes  ranged  from  4.67  (Golden  City  Lake)  to  7.24  (Grandview  Lake). 
Similarly  alkalinity  ranged  from  -1.1  (Golden  City)  to  9.48  mg-L"^  (Grandview  Lake;  see 
Table  2).  The  concentrations  of  DOC  ranged  from  2.7  mg-L^  in  Armishaw  Lake  to  10.5 
mg-L"^  in  Haller  Lake.  Values  of  TKN  ranged  from  200  (Armishaw)  to  566  A^g-L'^  (Moot 
Lake).  Total  phosphorus  ranged  from  4.7  (Grandview)  to  44.6  ^g-L"^  (Bittern  Lake). 
Aluminum  concentrations  ranged  from  0  /^g-L"^  (Bigwind)  to  128  /ug-L"^  (Axe  Lake). 
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Sulphate  values  ranged  from  3.95  (Fawn  2  Lake)  to  7.92  mg-L'^  (Chub  Lake). 

Watershed  2EC  -  Severn  River 

Lakes  in  this  watershed  empty  into  Georgian  Bay  and  Lake  Huron  through  the  Severn 
River.  We  sampled  6.7%  of  the  lakes  in  the  10  - 100  ha  size  range  (i.e.,  13  of  a  total  of 
193;  Cox  1978).  The  mean  area  of  these  lakes  was  34.8  ha  with  a  range  of  9.1  to  77 
ha  (Table  1).  Two  larger  lakes  were  also  sampled:  Clear  Lake  (103.8  ha)  and  Big  East 
Lake  (146.4  ha).  The  mean  depth  of  all  of  these  lakes  ranged  from  3.1  to  16.7  m. 
Elevation  ranged  from  318  to  404  mASL. 

In  this  watershed,  pH  ranged  from  5.11  (Nehemiah  Lake)  to  6.88  (Red  Chalk  Main), 
whereas  alkalinity  ranged  from  0.69  (Nehemiah  Lake)  to  5.14  mg-L'^  (Big  Orillia  Lake; 
see  Table  2).  The  DOC  ranged  from  2.1  (Clear  Lake)  to  9.5  mg-L""  (Nehemiah  Lake). 
Concentrations  of  TKN  ranged  from  180  (Red  Chalk  Main)  to  600  A^g-L''  (Nehemiah 
Lake).  Total  phosphorus  values  ranged  from  2.3  (Clear  Lake)  to  13.5  /ug-L"^  (Cinder 
Lake).  Aluminum  concentrations  ranged  from  4  (Clear  Lake)  to  200  /^g-L"^  (Nehemiah 
Lake).  Sulphate  values  ranged  from  6.3  (Big  Orillia  Lake)  to  8.7  mg-L"'  (Saw  Lake). 

Watershed  2HF  -  Cameron  Lake  Drainage 

Lakes  in  this  watershed  flow  into  Lake  Ontario  through  the  Cameron  River.  We 
sampled  20  lakes  in  this  watershed,  and  of  these,  only  7  were  between  10  and  99 
hectares  (mean  =  51.8,  range  =  16.0  to  98.7;  see  Table  1).  The  remaining  13  lakes 
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ranged  from  1 14.8  ha  to  4665.0  ha.  Most  of  these  larger  lakes  were  sampled  early  in 
the  study  before  the  final  lake-selection  criteria  were  established.  The  mean  depth  of 
these  lakes  ranged  from  4.7  to  31 .6  m  and  lake  elevation  ranged  from  307  to  41 1 
mASL. 

For  these  lakes,  the  pH  ranged  from  5.93  (Lipsy  Lake)  to  7.39  (Lower  Welch  Lake), 
whereas  alkalinity  ranged  from  0.63  (Goodwin  Lake)  to  40.94  mg-L"^  (Balsam  Lake;  see 
Table  2).  The  DOC  concentrations  ranged  from  2.0  (Kelly  Lake)  to  4.3  mg-L"^ 
(Basshaunt  Lake).  Total  phosphorus  values  ranged  from  2.0  (Boshkung  Lake)  to  13.0 
^ig-L'^  (Balsam).  Aluminum  concentrations  ranged  from  6  (Kelly  and  Margaret  lakes)  to 
54.0  /ugL'^  (Lipsy  Lake).    Sulphate  values  ranged  from  6.1  (Koshlong  Lake)  to  8.8 
mg-L'^  (Boshkung  Lake). 

Watersheds:  2HG  &  2KD  -  Scugog  River  &  Madawaska  River 

Because  of  the  small  number  of  lakes  that  were  sampled  in  these  two  watersheds,  we 
have  grouped  these  results  together.  Scugog  Lake  was  the  only  lake  sampled  in 
watershed  2HG  (Scugog  River).  This  lake  has  an  exceptionally  high  pH  (8.61)  and 
alkalinity  (118.9  mg-L'\  Table  2).  In  addition,  nutrient  concentrations  were  much  higher 
(i.e.,  TKN:  840  iug-L'\  TP:  46.2  /wg-L'^).  The  water  chemistry  of  Scugog  Lake  is  typical 
of  the  hardwater,  nutrient-enriched  lakes  of  the  Kawartha  Lakes  region  of  Ontario  (e.g., 
see  Berrill  1978).  We  plan  to  sample  additional  lakes  in  this  watershed  in  future  years. 

Two  lakes  were  sampled  in  watershed  2KD,  the  Upper  Madawaska  River  in  Algonquin 

20 


Provincial  Park  (i.e.,  Timberwolf  and  Louisa  lakes).  These  two  lakes  are  in  the  same 
watershed  as  Delano  and  Cradle  lakes  that  are  sampled  each  year  as  part  of  the 
temporal-lake  series  (see  David  et  al.  1994).  Although  Delano  and  Cradle  lakes  are 
relatively  small  (23.9  and  17.9  ha,  respectively),  Timberwolf  and  Louisa  are  larger  (167 
and  51 3  ha,  respectively;  Table  1 ).  The  water  chemistry  of  these  two  lakes  is  similar  to 
other  lakes  in  the  south-central  Ontario  watersheds  of  2EA,  2EB,  2EC  and  2HF  (Table 
2).  We  hope  to  sample  more  lakes  in  this  watershed  in  the  future. 

Fish  Community  Composition 

The  Ontario  Ministry  of  Natural  Resources  provided  inventory  data  on  the  fish  species 
present  in  many  of  the  lakes  in  this  survey  (Tables  3  and  4).  The  occurrence  of  crayfish 
predators,  such  as  smallmouth  bass  (Micropterus  dolomieui).  largemouth  bass  (M- 
salmoides).  rock  bass  (Ambloplites  rupestris).  yellow  perch  (Perca  flavescens).  walleye 
(Stizostedion  vitreum  vitreum).  brown  bullhead  (Ictalurus  nebulosus)  and  burbot  (Lota 
lota)  may  result  in  lower  crayfish  catches  in  baited  traps.  All  of  these  fish  species  could 
reduce  trap  catches  by  altering  crayfish  activity  patterns  and  related  behaviour  (Collins 
et  al.  1983,  Somers  and  Stechey  1986).  Alternatively,  these  species  and  others  such 
as  brook  trout  (Salvelinus  fontinalis)  and  lake  trout  (S.  namaycush)  may  directly  reduce 
crayfish  abundance  through  predation  (Momot  1967,  Stein  1977,  Somers  and  Green 
1993). 

As  above,  the  fish  species  inventories  are  organized  by  major  watershed.  The 
presence  of  smallmouth  bass,  largemouth  bass,  brown  bullhead  and  burbot  is  noted 
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because  these  species  are  probably  the  most  important  crayfish  predators  in  the  small 
lakes  of  south-central  Ontario  (Collins  et  al.  1983).  Other  fish  species,  such  as  yellow 
perch  and  rock  bass,  are  generally  too  small  to  be  major  crayfish  predators,  although 
this  is  not  the  case  in  the  shallower,  warmer  lakes  in  the  United  States  (Stein  1977, 
DiDonato  and  Lodge  1993). 

Watershed  2CF  -  Georgian  Bay  Island  Tributaries 

No  fish  species  were  collected  in  Clearwater  and  Lohi  lakes  (Table  4).  By  contrast. 
Middle  and  Hannah  lakes  have  recovered  in  the  last  decade  to  now  support  5  and  4 
fish  species,  respectively.  The  most  recent  surveys  indicate  that  Hannah  Lake  has 
creek  chub  (Semotilis  atromaculatus).  golden  shiner  (Notemigonus  cn/soleucas).  Iowa 
darter  (Etheostoma  exile)  and  yellow  perch,  whereas  golden  shiner,  Iowa  darter,  johnny 
darter  (E.  nigrum),  northern  redbelly  dace  (Phoxinus  eos)  and  yellow  perch  are  found  in 
Middle  Lake  (Gunn  and  Keller  1990,  Yan  et  al.  1996).  The  absence  offish  species  that 
are  important  crayfish  predators  suggests  that  the  crayfish  catches  in  these  lakes  are 
not  biased  by  the  influence  of  predators.  By  contrast,  most  lakes  in  the  other 
watersheds  support  a  variety  of  fish  species  that  are  important  crayfish  predators. 

Watershed  2EA  -  Georgian  Bay  Tributaries 

Fisheries  inventory  data  was  limited  to  the  major  species  of  sport  fish  in  several  of  the 
27  lakes  that  were  sampled  in  this  watershed.  Despite  these  limited  data,  only  Round, 
Back,  Frank  and  Grass  lakes  did  not  have  bass  or  bullhead  populations  (Table  4). 


22 


Smallmouth  bass  were  present  in  15  lakes,  largemouth  bass  were  collected  In  8  lakes, 
rock  bass  were  present  In  6  lakes,  and  brown  bullheads  were  found  In  13  lakes.  By 
contrast,  burbot  were  only  present  in  Clam  Lake.  Lake  trout  were  caught  in  Fairholme, 
Fowke,  Manitouwaba,  and  Grass  lakes.  These  4  lakes  all  have  maximum  depths 
greater  than  18  m  (see  Table  1)  supporting  the  observation  that  lake  trout  are  only 
found  in  relatively  deep  lakes  through  the  southern  part  of  their  range  in  Canada  (Scott 
and  Crossman  1973,  Gunn  et  al.  1988).  The  occurrence  of  crayfish  predators  (i.e., 
bass  or  bullheads)  in  85%  of  the  surveyed  lakes  in  this  watershed  suggests  that 
crayfish  CUE  may  be  reduced  relative  to  lakes  in  other  watersheds. 

Watershed  2EB  -  Moon  River  and  Go  Home  River 

Survey  data  were  available  for  all  but  two  of  the  29  lakes  that  were  sampled  in  this 
watershed  (i.e..  Groves  and  Haller  lakes;  Table  4).  Smallmouth  bass  were  found  in  20 
of  these  27  lakes.  Largemouth  bass  were  present  in  12  lakes  and  bullheads  were 
present  in  17  lakes.  Burbot  were  only  caught  in  Grandview  Lake.  Although  none  of 
these  crayfish  predators  was  present  in  Armishaw,  Clayton,  Langford  or  Lynx  lakes, 
crayfish  predators  were  found  in  86%  of  the  29  lakes.  As  a  result,  crayfish  trap  catches 
may  be  somewhat  lower  than  catches  from  lakes  in  those  watersheds  where  crayfish 
predators  are  less  common. 

Watershed  2EC  -  Severn  River 

Survey  data  were  available  for  all  15  lakes  that  were  sampled  in  this  watershed. 
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Bullheads  were  present  in  6  of  these  15  lakes  (Table  4).  Smallmouth  and  largemouth 
bass  were  found  in  8  lakes.  Rock  bass  were  only  caught  in  Cinder  Lake.  Major 
crayfish  predators  were  absent  from  Cinder,  Nehemiah,  Saw  and  Shoe  lakes.  Burbot 
or  lake  trout  were  only  present  in  Red  Chalk  and  Clear  lakes.  Consequently,  73%  of 
the  surveyed  lakes  in  this  watershed  had  major  crayfish  predators.  This  proportion  is 
markedly  lower  than  the  85-90%  values  for  the  2EA,  2EB  and  2HF  watersheds,  and 
hence  crayfish  catches  may  be  somewhat  higher  than  in  the  other  watersheds. 

Watershed  2HF  -  Cameron  Lake  Drainage 

Fish  survey  data  were  available  for  all  20  lakes  that  were  sampled  in  the  2HF 
watershed  (Table  4).  Only  4  lakes  (i.e.,  Kelly,  Lipsy,  Little  Kennisis  and  Margaret)  did 
not  have  smallmouth  bass  populations.  Of  these  4  lakes,  Little  Kennisis  supported 
bullhead  and  rock  bass  populations.  A  total  of  10  lakes  had  bullhead  populations.  In 
contrast  with  the  other  watersheds,  rock  bass  were  more  common  in  these  lakes.  That 
is,  rock  bass  were  caught  in  10  of  the  20  lakes.  Largemouth  bass  were  found  in  7  lakes 
and  burbot  were  caught  in  Balsam,  Big  Hawk,  Kushog,  Sherborne,  St.  Nora  and  Twelve 
Mile  lakes.  Lake  trout  were  present  in  16  of  the  lakes.  In  total,  18  lakes  (i.e.,  90%) 
supported  populations  of  bullheads  or  bass.  As  a  result,  crayfish  catches  from  these 
lakes  may  be  somewhat  lower  because  of  the  prevalence  of  predatory  fish. 

Watershed  2HG  -  Scugog  River 

Scugog  Lake  was  the  only  lake  that  was  sampled  in  the  Scugog  River  watershed  (Table 
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4).  This  lake  supports  populations  of  several  predatory  fish  species  including: 
muskellunge  (Esox  masquinongy).  smallmouth  bass,  largemouth  bass,  rock  bass, 
yellow  (Ictalurus  natalis)  and  brown  bullheads,  and  walleye.  Many  other  fish  species 
also  occur  in  this  lake,  but  most  of  these  species  are  not  believed  to  be  significant 
crayfish  predators. 

Watershed  2KD  -  Upper  Madawaska  River 

Only  two  lakes  in  the  2KD  watershed  were  sampled  during  this  study.  Fisheries  data 
was  available  for  both  of  these  lakes  (Table  4).  Timberwolf  Lake  contained  populations 
of  brook  trout,  lake  trout,  burbot,  brown  bullhead  and  yellow  perch,  as  well  as  a  number 
of  other  fish  species  that  likely  have  little  impact  on  crayfish.  By  contrast,  Louisa  Lake 
supported  brook  and  lake  trout,  but  not  brown  bullhead  or  bass  populations.  Although 
some  crayfish  predators  were  captured,  the  impact  offish  predation  on  crayfish 
populations  in  these  two  lakes  is  probably  less  than  in  lakes  in  the  other  watersheds. 

Crayfish  Catch  Data 

Nine  species  of  crayfish  occur  in  Ontario  (Crocker  and  Barr  1968).  The  5  native 
species  that  were  collected  in  this  study  probably  followed  the  receding  glaciers 
northwards  (i.e.,  O.  propinquus.  O.  virilis.  O.  immunis,  C.  bartoni  and  C.  robustus: 
Crocker  and  Barr  1968,  Guiasu  et  al.  1996).  By  contrast,  2  species,  O.  rusticus  and  O. 
obscurus.  were  probably  introduced  into  Ontario  waters  by  fishermen  (Crocker  and  Barr 
1968,  Berrill  1978,  Corey  1988,  Momot  1992). 
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Quality  Assurance/Quality  Control 

Seven  lakes  were  sampled  twice  during  this  study  to  estimate  the  degree  of  year-to- 
year  variation  in  crayfish  relative  abundance  (Table  5).  The  length  of  time  between 
sampling  visits  varied  from  3  to  5  years.  Bigwind,  Dickie,  Chub  and  Moot  lakes  were 
sampled  in  1991,  whereas  Twelve  Mile  Lake  was  surveyed  in  1990  and  Red  Chalk 
Lake  (with  separate  East  and  Main  basins)  was  originally  surveyed  in  1989.  All  of  these 
lakes  were  re-surveyed  in  1994. 

Year-to-year  variation  in  the  Quality  Assurance/Quality  Control  (QA/QC)  data  was 
evident  (Table  5),  although  no  consistent  trends  were  revealed  (Figure  2).  The  catches 
of  C.  bartoni.  O.  propinquus  and  O.  virilis  were  scattered  on  both  sides  of  the  1 :1  line 
suggesting  that  none  of  the  catches  of  any  particular  species  was  biased  in  a  consistent 
way.  The  hypothesized  1:1  relationship  with  an  intercept  through  the  origin  (i.e.,  0,  0) 
accounted  for  37.5%  of  the  variation  in  the  CUE  for  the  second  visit  (i.e.,  r^=0.375). 
Undoubtedly  the  large  catch  of  O.  propinquus  from  Bigwind  Lake  in  1994  (i.e.,  8.68 
versus  3.80  in  1991;  Table  5)  contributed  a  sizeable  amount  to  the  unexplained 
variation  in  the  1994  CUE. 

Three  lakes  had  relatively  large  trap  catches  on  each  sampling  date  (i.e.,  Bigwind,  Red 
Chalk  East  and  Red  Chalk  Main;  Table  5,  Figure  2).  In  Bigwind  Lake,  the  CUE  of  O. 
propinquus  and  O.  virilis  increased  from  1991  to  1994,  whereas  the  catch  of  C.  bartoni 
fell  to  zero  in  1994.  This  zero  catch  does  not  indicate  that  C.  bartoni  is  now  extinct  in 
Bigwind  Lake  (cf.  France  and  Collins  1993)  because  we  are  unable  to  distinguish  a  true 
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extinction  from  a  relative  abundance  that  is  below  the  detection  limit  for  this  sampling 
protocol.  Until  appropriate  detection  limits  are  established,  we  recommend  caution 
when  interpreting  crayfish  CUE  values  that  approach  zero.  That  is,  the  disappearance 
and  subsequent  re-appearance  of  particular  taxa  may  simply  reflect  relative 
abundances  that  are  near  or  below  the  detection  limits.  The  time-series  or  temporal 
study  by  David  et  al.  (1994)  may  also  reflect  this  problem  because  the  crayfish  catches 
from  several  lakes  also  dropped  to  zero. 

All  of  the  crayfish  species  caught  in  Red  Chalk  Lake  in  1989  were  also  caught  in  1994 
(Table  5).  However,  the  catches  of  O.  propinquus  and  O.  virilis  declined  in  both  Red 
Chalk  East  and  Red  Chalk  Main.  The  catches  of  C.  bartoni  in  Red  Chalk  Main 
increased.  No  cambarids  were  caught  in  Red  Chalk  East,  despite  their  presence  in  the 
main  basin  of  Red  Chalk  Lake. 

The  remaining  QA/QC  lakes  had  very  low  catches  (Table  5).  For  example,  the  catch  of 
C.  bartoni  in  Moot  Lake  changed  very  little  from  1991  (0.13)  to  1994  (0.11).  The  CUE 
of  O.  virilis  in  Chub  Lake  was  quite  low  in  1991  (i.e.,  0.04)  and  fell  to  zero  in  1994,  much 
like  the  catch  of  C.  bartoni  in  Bigwind  Lake.  In  addition,  no  crayfish  were  caught  in 
Twelve  Mile  Lake  in  1990,  although  O.  virilis  were  captured  by  SCUBA  divers  during 
that  same  period  (M.  Berrill  and  G.  Taylor,  Trent  University,  unpubl.  data).  This 
discrepancy  supports  the  hypothesis  that  the  presence  of  predatory  fish  alters  crayfish 
trappability  (Collins  et  al.  1983,  Somers  and  Green  1993)  because  Twelve  Mile  Lake 
supports  populations  of  rock  bass,  burbot,  smallmouth  and  largemouth  bass. 
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Subsequent  re-sampling  of  Twelve  Mile  Lake  with  baited  traps  in  1994  resulted  in  the 
collection  of  several  O.  propinquus  and  O.  virilis.  By  contrast,  no  crayfish  were  caught 
in  Dickie  Lake  on  either  visit. 

The  QA/QC  results  emphasize  the  need  for  temporal  (or  time-series)  data  in  order  to 
observe  and  quantify  natural  fluctuations  in  crayfish  abundances  (e.g.,  David  et  al. 
1994).  More  importantly,  we  should  establish  detection  limits  to  avoid  reporting  the 
absence  of  crayfish  in  lakes  where  the  populations  are  present,  but  at  low  relative 
abundances.  Although  comparisons  spanning  several  years  are  useful,  we  also  hope 
to  re-sample  several  lakes  to  quantify  both  seasonality  (e.g.,  Somers  and  Green  1993) 
as  well  as  the  repeatability  of  our  CUE  estimates  (e.g.,  Reid  et  al.  1995). 

Male  -  Female  Trapping  Biases 

Baited  traps  catch  proportionately  more  large  male  crayfish  (Momot  et  al.  1978,  Somers 
and  Stechey  1986).  Because  of  this  recognised  bias  in  trap  data,  the  catch  of  females 
is  often  ignored  when  estimating  crayfish  CUE.  This  approach  usually  results  In  simply 
reporting  the  mean  number  of  adult  males  per  trap  (Capelli  and  Magnuson  1983, 
France  1993).  However,  Malley  and  Reynolds  (1979)  have  also  cautioned  that  trap 
catches  provide  biased  estimates  of  population  structure,  including  size  distributions 
and  sex  ratios.  Although  we  recognise  these  concerns  (e.g.,  Somers  and  Stechey 
1986,  Somers  and  Green  1993),  we  report  and  analyse  CUE  data  for  both  males  and 
females  combined  because  of  the  relatively  small  catches  in  most  lakes  (e.g.,  see 
Table  5),  and  because  crayfish  were  only  collected  during  the  mid-summer  months 
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when  trap  catches  are  less  biased  (Somers  and  Green  1993). 

Year-to-year  differences  in  the  male-to-female  catch  ratios  were  evident  in  this  survey, 
although  no  consistent  patterns  were  apparent  (Table  5).  For  example,  more  male  O. 
virilis  were  caught  in  1989,  1991  and  1994,  whereas  more  females  were  caught  in  1990 
and  1992  (Table  6).  Because  recently  moulted  crayfish  are  rarely  caught  in  baited 
traps,  these  differences  are  probably  associated  with  annual  differences  in  crayfish 
moulting  patterns  that  are  influenced  by  water  temperature  (Somers  and  Green  1993). 
Males  and  females,  as  well  as  the  different  species,  moult  at  slightly  different  times  in  a 
given  year  (Crocker  and  Barr  1968,  Corey  1988,  Somers  and  Green  1993).  As  a  result, 
it  is  not  surprising  that  the  sex  ratios  differ  among  taxa  in  different  years.  Because  the 
seasonal  pattern  in  water  temperature  probably  affects  the  proportion  of  males  in  the 
total  catch,  long-term  changes  in  sex  ratios  of  crayfish  catches  may  reflect  changes 
associated  with  global  warming.  This  relationship  should  be  examined  as  additional 
data  are  added  to  the  temporal  data  series. 

In  the  set  of  temporal  reference  lakes,  David  et  al.  (1994)  reported  more  males  than 
females  for  trap  catches  of  C.  bartoni  and  O.  virilis.  but  not  O.  propinquus.  We  found 
similar  results  from  this  set  of  100  lakes  (Tables  5  and  6).  The  larger  catches  of  males 
in  many  lakes  may  reflect  the  presence  of  predatory  fish  (Collins  et  al.  1983,  Somers 
and  Stechey  1986,  Somers  and  Green  1993).  For  example,  Stein  (1977)  found  that 
male  crayfish  were  less  vulnerable  to  fish  predation  than  females.  As  a  result,  male 
crayfish  tend  to  be  more  active  than  females  in  lakes  with  predators  (Collins  et  al.  1983, 
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Somers  and  Stechey  1986).  Thus,  males  are  more  likely  to  encounter  traps  in  lakes 
with  predatory  fish,  and  consequently  the  sex  ratios  of  crayfish  catches  from  these 
lakes  may  be  biased  towards  larger  numbers  of  males.  The  absence  of  a  male-biased 
sex  ratio  in  the  catches  of  O.  propinquus  may  arise  because  this  species  is  generally 
smaller  in  overall  size,  and  consequently,  both  sexes  may  be  equally  susceptible  to 
predation  by  fish. 

Size  Distributions  in  the  Trap  Catches 

Cambarus  robustus  was  the  largest  crayfish  caught  during  this  study,  with  mean 
carapace  lengths  ranging  from  35  to  48  mm  (Table  7).  By  contrast,  the  smallest  mean 
carapace  lengths  were  obtained  for  O.  propinquus.  which  ranged  from  18  to  35  mm. 
These  results  are  not  surprising  since  C.  robustus  is  the  largest  crayfish  and  O. 
propinquus  is  the  smallest  crayfish  found  in  Ontario  (Crocker  and  Barr  1968). 

Ranges  in  carapace  lengths  for  the  catches  of  each  species,  as  well  as  overall  means 
for  males  and  females  were  tabulated  by  watershed  (see  Table  7).  In  general,  the 
mean  values  for  a  given  species  were  similar  across  watersheds.  The  exceptions 
tended  to  be  those  watersheds  where  relatively  few  lakes  were  sampled.  As  in  the 
temporal-lakes  study  (David  et  ai.  1994),  Form  I  males  were  usually  the  largest  crayfish 
caught,  whereas  females  and  Form  11  males  were  slightly  smaller  on  average. 

France  (1985)  suggested  that  crayfish  growth  rates,  and  hence  size  distributions,  were 
a  function  of  lake  productivity.  Since  carapace  lengths  were  recorded  for  most  of  the 
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crayfish  that  were  caught,  we  originally  tabulated  the  minimum,  maximum  and  mean 
carapace  lengths  for  each  species  in  each  lake.  However,  the  relatively  small  sample 
sizes  from  many  lakes  limited  the  value  of  these  comparisons.  After  subdividing  the 
catches  into  females,  reproductive  males  and  non-reproductive  males  for  each  species 
in  each  lake  we  often  found  fewer  than  10  animals  in  many  categories.  As  a  result, 
minimum,  mean  and  maximum  sizes  are  probably  not  representative  for  many  catches. 

Relatively  few  of  the  lakes  yielded  catches  with  more  than  10  individuals  of  a  particular 
species  or  sex.  For  example,  only  the  catches  of  C.  bartoni  from  10,  6  and  4  lakes 
resulted  in  10-or-more  females.  Form  I  males  or  Form  II  males,  respectively.  The  data 
for  O.  propinquus  were  similar,  with  10-or-more  females,  reproductive  and  non- 
reproductive  males  from  8,  4  and  6  lakes,  respectively.  On  average,  larger  numbers  of 
O.  virilis  were  captured  in  the  various  lakes.  Ten-or-more  individuals  were  collected 
from  16,  10  and  14  lakes,  respectively.  The  largest  Pearson  correlations  between 
average  female  carapace  length  and  the  various  water-chemistry  parameters  involved 
K  (0.759),  TKN  (0.690)  and  Al  (0.673).  The  largest  correlations  for  the  Form  I  males 
were  with  K  (0.708),  Mg  (0.678)  and  alkalinity  (0.661).  By  contrast,  correlations  for  the 
Form  II  males  were  lower,  although  the  highest  3  values  involved  the  same  variables  as 
the  female  data  (i.e.,  K:  0.500,  Al:  0.497,  and  TKN:  0.474).  Given  the  limited  number  of 
lakes,  we  caution  against  over-interpreting  these  results.  However,  the  surprising 
consistency  of  the  correlations  suggests  that  further  work  in  this  area  may  prove  to  be 
rewarding. 
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Crayfish  Relative  Abundances 

Like  other  animals,  crayfish  populations  are  subject  to  fluctuations  in  abundance  from 
one  year  to  the  next  (Lodge  et  al.  1986,  David  et  al.  1994).  Because  absolute 
abundance  extirpates  based  on  crayfish  mark-recapture  studies  are  exceedingly  rare 
(e.g.,  Emery  1975),  trap  catches  are  frequently  used  to  provide  estimates  of  crayfish 
relative  abundance  (e.g.,  Capelli  and  Magnuson  1983,  Lodge  et  al.  1986).  However, 
trap  catches  are  affected  by  many  factors  including  the  presence  of  predatory  fish,  bait 
type  and  seasonal  cues  that  influence  crayfish  life  history  (e.g.,  Collins  et  al.  1983, 
Somers  and  Stechey  1986,  Somers  and  Green  1993).  As  a  result,  the  presence  or 
absence  of  a  particular  crayfish  species  in  a  survey  can  only  be  regarded  as  a 
"snapshot"  estimate  for  that  sampling  period.  Quality  Assurance/Quality  Control 
procedures  (e.g.,  Reid  et  al.  1995)  help  to  assess  the  accuracy  and  precision  of  the 
data,  and  thereby  provide  a  measure  of  confidence  for  relative  abundance  estimates. 
In  this  context,  time-trend  data  (e.g.,  David  et  al.  1994)  and  replicated  sampling  provide 
initial  estimates  of  CUE  reliability  and  repeatability  (see  Figure  2). 

The  time-trend  data  from  14  south-central  Ontario  lakes  revealed  that  crayfish 
abundances  in  several  lakes  declined  from  1988  to  1994  (David  et  al.  1994).  In  some 
cases,  species  composition  also  changed  from  year  to  year  when  particular  species 
were  not  caught.  In  other  lakes,  crayfish  abundances  changed  very  little.  Lodge  et  al. 
(1986)  present  evidence  of  pronounced  temporal  changes  in  crayfish  abundances  after 
O.  rusticus  was  introduced  into  a  Wisconsin  lake.  Because  of  these  types  of 
fluctuations  in  relative  abundance,  we  propose  that  this  spatial  survey  should  be 
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repeated  in  order  to  accurately  model  the  responses  of  crayfish  populations  to 
anthropogenic  stress.  Further  sampling  is  also  necessary  in  order  to  distinguish  the 
impacts  of  introduced  species  (i.e.,  biological  perturbations)  from  other  man-made 
disturbances.  This  distinction  is  relevant  because  there  are  no  acceptable  remedial 
measures  to  reverse  the  impacts  of  introduced  species  (Momot  1992).  By  contrast, 
there  are  numerous  options  to  reverse  physical  and  chemical  alterations  (e.g.,  see 
Gunn  and  Keller  1990,  Yan  et  al.  1996). 

Watershed  2CF  -  Georgian  Bay  Island  Tributaries 

Only  one  crayfish  species  was  caught  in  the  4  lakes  surveyed  in  this  watershed  (Table 
8).  Cambarus  bartoni  was  found  in  Lohi  and  Clearwater  lakes,  but  no  crayfish  were 
caught  in  Middle  and  Hannah  lakes.  Despite  extreme  chemical  conditions  (Sudbury 
Environmental  Study  1982),  crayfish  populations  inhabited  CleanA/ater  and  Lohi  lakes  in 
the  early  1970s  when  the  pH  was  4.2.  There  is  no  comparable  information  for  Middle 
and  Hannah  lakes.  Cambarus  bartoni  is  an  acid-tolerant  species  (Berrill  1978,  Berrill  et 
al.  1985),  so  its  occurrence  in  low-pH  waters  is  not  surprising.  Whether  or  not  C. 
robustus.  O.  virilis  or  O.  propinquus  once  inhabited  this  watershed  is  unknown,  although 
attempts  to  re-introduce  crayfish  into  several  Sudbury-area  lakes  are  ongoing  (John 
Gunn,  MNR  Sudbury,  pers.  comm.). 

Watershed  2EA  -  Georgian  Bay  Tributaries 

Six  crayfish  species  were  found  in  the  Georgian  Bay  Tributaries  watershed  (Table  8). 
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In  those  lakes  where  crayfish  were  captured,  the  CUE  (catch  per  unit  effort)  ranged 
from  0.02  to  1 .63.  The  average  catch  for  all  27  lakes  was  0.362  crayfish  per  trap  per 
night  (see  Tables  5  and  8).  Cheer  Lake  was  the  only  lake  where  no  crayfish  were 
caught  (i.e.,  3.7%  of  the  lakes).  This  lake  is  relatively  large  (area  =  90.1  ha)  so  the 
absence  of  crayfish  Is  somewhat  surprising,  especially  since  the  water  chemistry  is  not 
unusual  with  respect  to  other  lakes  in  the  watershed  (i.e.,  pH  =  5.92,  A!  =  39  mQ-L '^  see 
Table  2).  In  all  likelihood  crayfish  may  be  present  in  this  lake,  but  their  density  must  be 
low  relative  to  the  other  lakes. 

Orconectes  propinquus  was  the  most  common  crayfish  species  caught  in  this 
watershed,  occurring  in  1 7  of  27,  or  64%  of  the  lakes  (Table  8).  In  two  of  these  lakes 
(Bat  and  Richmond),  O.  propinquus  was  the  only  crayfish  species  caught.  By  contrast, 
O.  immunis  was  the  rarest  species,  occurring  in  only  4  lakes  (Old  Mans,  Limestone, 
Fowke  and  Frank).  Surprisingly,  O.  immunis  was  found  in  Frank  Lake,  a  relatively  small 
(19.4  ha)  high-elevation  lake  (389  mASL,  see  Table  1).  By  contrast,  the  other  3  lakes 
were  relatively  large  lakes  at  elevations  below  300  mASL.  Just  how  O.  immunis 
colonized  Frank  Lake,  but  not  the  lakes  immediately  downstream  remains  a  mystery. 
Perhaps  this  population  was  introduced  by  fishermen. 

The  introduced  species  O.  obscurus  was  found  in  6  lakes  in  this  watershed  (Bells, 
Crawford,  Old  Mans,  Simmons,  Black  and  Himbury;  Table  8).  This  species  was 
probably  introduced  into  this  watershed  by  fishennen,  since  the  normal  geographic 
range  of  this  crayfish  is  eastern  Ohio,  Pennsylvania  and  New  York  (Crocker  and  Barr 
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1968).  The  presence  of  sportfish  in  these  lakes  lends  credence  to  this  hypothesis  (see 
Table  4). 

Cambarus  bartoni  and  C.  robustus  were  caught  in  6  and  7  lakes,  respectively  (Table  8). 
Cambarus  bartoni  was  relatively  rare,  ranging  in  abundance  from  0.02  to  a  maximum  of 
0.52  animals  per  trap  night.  By  contrast,  C.  robustus  was  quite  common  in  Fowke  Lake 
with  a  CUE  of  1 .48.  In  the  remaining  6  lakes,  the  abundance  of  C.  robustus  ranged 
from  0.02  to  0.43,  much  like  C.  bartoni.  Orconectes  virilis  was  also  common  in  these 
lakes  with  a  CUE  ranging  from  0.02  to  0.41  based  on  12  lakes  (44%). 

Watershed  2EB  -  Moon  River  and  Go  Home  River 

Four  crayfish  species  were  found  in  the  31  lakes  that  were  surveyed  in  this  watershed: 
C.  bartoni.  C.  robustus.  O.  propinquus  and  O.  virilis  (Table  8).   No  introduced  species 
was  captured.  The  most  common  species,  O.  virilis,  was  caught  in  13  or  42%  of  the 
lakes.  Orconectes  propinquus  was  collected  in  12  (39%)  of  the  lakes,  whereas  C. 
robustus  was  the  rarest  species,  occurring  in  only  two  lakes.  By  contrast,  C.  bartoni 
was  found  in  8  or  26%  of  the  lakes. 

Although  no  crayfish  were  caught  in  42%  (i.e.,  13)  of  the  lakes  in  this  watershed,  the 
average  CUE  was  almost  twice  the  average  value  for  watershed  2EA  (i.e.,  0.695;  see 
Table  8).  In  those  lakes  where  crayfish  were  found,  the  total  CUE  ranged  from  0.02  to 
12.78.  The  CUE  of  12.78  was  the  largest  catch  in  any  of  the  lakes,  representing  a  total 
of  690  C.  bartoni  in  54  traps  in  one  night.  The  two  lakes  with  the  highest  CUE  (i.e., 
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Bigwind  [9.94]  and  Clayton  [12.78])  do  not  have  bass,  bullhead  or  burbot  populations. 
As  noted  earlier,  the  absence  of  predatory  fish  may  contribute  to  these  large  catches. 

Despite  the  large  numbers  of  C.  bartoni  caught  in  Clayton  Lake,  the  catches  of  C. 
bartoni  in  the  other  lakes  were  quite  small  ranging  from  0.02  to  0.13  (Table  8). 
Similarly,  the  catches  of  C.  robustus  were  only  0.08  and  0.18  crayfish  per  trap  night.  By 
contrast,  the  catches  of  O.  propinquus  ranged  from  0.02  to  8.68,  whereas  the  catches 
of  O.  virilis  ranged  from  0.02  to  1.26.  Based  on  these  results,  the  orconectids  had 
higher  relative  abundances  in  the  lakes  in  the  2EB  watershed  relative  to  the  2EA 
watershed;  however,  far  more  of  the  2EB  lakes  were  devoid  of  crayfish  (i.e.,  42% 
versus  4%). 

in  watershed  2EA,  C.  bartoni  was  caught  in  22%  of  the  lakes,  whereas  O.  propinquus 
and  O.  virilis  were  caught  in  64%  and  44%  of  the  lakes,  respectively.  Cambarus  bartoni 
was  found  in  26%  of  the  lakes  in  the  2EB  watershed.  By  contrast,  O.  propinquus  was 
caught  in  39%,  whereas  O.  virilis  was  collected  in  42%  of  the  2EB  lakes.  Despite  the 
relatively  large  proportion  of  2EB  lakes  without  crayfish,  the  frequencies  of  C.  bartoni 
and  O.  virilis  were  similar  to  values  for  the  2EA  lakes.  Only  O.  propinquus  seemed  to 
be  under-represented  in  the  2EB  watershed.  Berrill  et  al.  (1985)  found  that  O. 
propinquus  is  intolerant  of  conditions  characteristic  of  acidic  lakes.  As  a  result, 
populations  of  O.  propinquus  may  have  been  lost  from  lakes  in  the  2EB  watershed 
much  like  the  population  of  O.  virilis  in  Plastic  Lake  in  watershed  2HF  (e.g.,  France  and 
Collins  1993,  David  et  al.  1994). 
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Simple  comparisons  between  lakes  in  the  2EA  and  2EB  watersheds  do  not  justify  the 
observed  differences  in  crayfish  catches  (Table  8).  Both  sets  of  lakes  span  similar 
ranges  in  lake  morphology  (Table  1)  and  water  chemistry  (Table  2).  However,  70%  of 
the  lakes  that  were  sampled  in  the  2EA  watershed  are  below  an  elevation  of  300 
mASL.  By  contrast,  only  23%  of  the  lakes  that  were  sampled  in  the  2EB  watershed  are 
below  300  mASL.  Because  high-elevation  lakes  tend  to  have  less  glacial  till  in  their 
watersheds,  the  buffering  capacity  of  these  lakes  also  tends  to  be  lower  (see  Table  2 
and  compare  alkalinity  and  conductivity;  Dillon  et  al.  1984).  Perhaps  historically,  the 
higher-elevation  lakes  in  the  2EB  watershed  were  affected  by  lake  acidification  to  a 
greater  extent  than  the  survey  lakes  in  the  2EA  watershed.  Physical  barriers  to  re- 
colonization  (e.g.,  the  elevation  gradient)  may  have  prolonged  the  time  required  for  the 
biological  recovery  of  these  2EB  lakes  (e.g.,  see  Gunn  et  al.  1988,  Gunn  and  Keller 
1990,  France  and  Collins  1993,  Yan  et  al.  1996).  In  addition,  the  formation  of  ice  on 
the  bottom  of  the  inter-connecting  streams  during  the  winter  may  limit  the  dispersal 
ability  of  crayfish  in  this  part  of  Ontario  (e.g.,  see  Aiken  1968). 

Alternatively,  the  post-glacial  dispersal  of  O.  propinquus  may  have  occurred  somewhat 
later  than  other  species  such  as  C.  bartoni  and  O.  virilis,  thereby  accounting  for  the 
reduced  frequency  of  occurrence  of  O.  propinquus  in  the  higher  elevation  lakes  in 
south-central  Ontario  (see  Crocker  and  Barr  1968).  However,  the  fact  that  O. 
propinquus  is  often  found  in  the  larger  lakes  above  an  elevation  of  300  mASL  (i.e., 
compare  Tables  1  and  8)  suggests  that  historical  acidification,  rather  than  post-glacial 
dispersal  is  implicated  (see  Figure  3). 
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Watershed  2EC  -  Severn  River 

Only  3  species  of  crayfish  were  caught  in  the  15  lakes  that  were  sampled  in  the  Severn 
River  watershed  (Table  8).  All  15  of  these  lakes  were  at  elevations  greater  than  300 
mASL.  As  a  result,  finding  relatively  few  crayfish  species  in  this  watershed  is  consistent 
with  the  hypothesis  that  elevation  affects  crayfish  species  distributions  in  south-central 
Ontario  (Crocker  and  Barr  1968;  and  see  Figure  3). 

All  3  species  (i.e.,  C.  bartoni.  O.  propinquus  and  O.  virilis)  co-occurred  in  5  lakes  (Table 
8).  No  crayfish  were  caught  in  3  (20%)  of  the  lakes,  whereas  10  lakes  contained  more 
than  one  crayfish  species.  Orconectes  virilis  was  the  most  common  crayfish,  occurring 
in  12  or  80%  of  the  lakes.  Cambarus  bartoni  was  found  in  8  (53%)  of  the  lakes, 
whereas  O.  propinquus  occurred  in  6  lakes  (i.e.,  40%).    By  comparison  to  the  2EA  and 
2EB  watersheds,  both  C.  bartoni  and  O.  virilis  were  found  in  twice  as  many  lakes  (as  a 
percentage)  in  the  2EC  watershed.  Orconectes  propinquus  was  caught  in  40  and  39% 
of  the  lakes  in  the  2EC  and  2EB  watersheds,  respectively.  These  values  are  markedly 
lower  than  their  occurrence  in  64%  of  the  lakes  that  were  sampled  in  the  2EA 
watershed. 

The  CUE  of  C.  bartoni  ranged  from  0.02  to  0.44,  much  like  values  for  the  lakes  in  the 
2EA  watershed,  but  higher  than  most  of  the  values  from  the  2EB  watershed  (Table  8). 
On  average,  the  total  CUE  for  the  2EC  lakes  was  1 .869  which  is  much  higher  than  the 
other  watersheds.  The  O.  propinquus  CUE  ranged  from  0.07  to  2.80.  These  values 
are  similar  to  the  catches  from  lakes  in  the  2EB  watershed,  but  they  are  higher  than  the 
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values  for  the  2EA  watershed.  Similarly,  the  CUE  of  O.  virilis  ranged  from  0.02  to  4.35 
animals  per  trap  night.  These  catches  are  noticeably  higher  than  the  CUE  of  O.  virilis  in 
lakes  in  the  2EB  and  2EA  watersheds. 

Simple  comparisons  between  the  lakes  in  the  Severn  River  watershed  (2EC)  and  the 
lakes  in  watersheds  2EA  and  2EB  emphasize  that  all  of  the  2EC  lakes  are  above  300 
mASL  (Table  1).  Lake  sizes  and  average  depths  span  similar  ranges  much  like  the 
water  chemistry  parameters  (Table  2).  One  surphsing  difference  is  the  fact  that  crayfish 
predators  represented  by  bullheads  and  smallmouth  and  largemouth  bass  are  found  in 
1 1  or  73%  of  the  1 5  lakes.  These  same  predators  are  found  in  85  and  86%  of  the  2EA 
and  2EB  lakes,  respectively  (Table  4).  Perhaps  the  reduced  frequency  of  crayfish 
predators  contributes  to  the  larger  crayfish  catches  in  the  2EC  lakes.  Alternatively, 
some  unmeasured  parameter  may  be  involved. 

Watershed  2HF  -  Cameron  Lake  Drainage 

Six  crayfish  species  were  found  in  the  Cameron  Lake  watershed,  although  two  of  these 
species  (O.  immunis  and  O.  rusticus)  were  only  found  in  one  lake  (Table  8).  No 
crayfish  were  caught  in  9  lakes  (45%).  This  relatively  high  percentage  of  lakes  without 
crayfish  is  comparable  to  the  2EB  watershed  (42%),  although  this  percentage  is  much 
higher  than  the  values  for  the  2EA  (4%)  and  2EC  (20%)  watersheds. 

Six  of  the  20  lakes  (30%)  had  C.  bartoni  populations,  whereas  C.  robustus  and  O. 
propinquus  were  each  caught  in  3  lakes  (Table  8).  The  proportion  of  the  2HF  lakes  with 
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C.  bartoni  is  similar  to  the  2EA  (22%)  and  2EB  (26%)  watersheds,  but  lower  than  2EC 
(53%).  Similarly,  the  relative  frequency  of  O.  virilis  (8  lakes  or  40%)  is  consistent  with 
the  2EA  (44%)  and  2EB  (42%)  watersheds,  but  not  2EC  (80%).  By  contrast,  O. 
propinquus  was  caught  in  63%  of  the  2EA  lakes,  40%  of  the  2EB  and  2EC  lakes,  but 
only  15%  of  the  2HF  lakes.  Some  of  these  differences  may  be  attributed  to  the  fact  that 
all  but  one  of  the  2HF  lakes  were  at  elevations  above  300  mASL. 

In  those  lakes  where  crayfish  were  caught,  the  total  CUE  ranged  from  0.01  to  3.36 
crayfish  per  trap  per  night  (Table  8).  The  average  CUE  was  0.409  which  is  similar  to 
the  2EA  watershed  (0.362),  but  lower  than  the  2EB  (0.695)  and  2EC  (1 .869) 
watersheds.  The  catches  of  C.  bartoni  were  generally  quite  low,  ranging  from  0.01  to 
0.22,  with  a  single  high  value  of  1.09  in  Little  Buckhorn  Lake.  The  3  catches  of  O. 
propinquus  ranged  from  0.02  to  0.39.  By  contrast,  the  CUE  of  O.  virilis  was  somewhat 
higher,  ranging  from  0.04  to  a  maximum  of  2.21  in  Basshaunt  Lake. 

Eight  of  the  lakes  that  were  surveyed  in  the  Cameron  Lake  watershed  had  more  than 
one  species  of  crayfish  (Table  8).  Balsam  Lake  was  the  only  lake  where  O.  rusticus 
was  found.  Orconectes  propinquus.  O.  virilis  and  C.  robustus  were  also  collected  in 
Balsam  Lake.  As  previously  mentioned,  O.  rusticus  has  been  introduced  into  Ontario 
from  the  United  States.  Crocker  and  Barr  (1968)  originally  found  it  in  6  locations  in 
Ontario.  Since  that  time  O.  rusticus  has  been  reported  from  the  Kawartha  Lakes  area 
(Berrill  1978)  and  near  Thunder  Bay  (Momot  1992).  Berrill  (1978)  proposed  that  O. 
rusticus  may  compete  with  O.  propinquus  and  O.  virilis.  This  hypothesis  has  been 
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substantiated  by  field  studies  where  O.  rusticus  has  displaced  O.  propinquus  and  O. 
virilis  in  several  lakes  (Capelli  1982,  Lodge  et  al.  1986,  Olsen  et  al.  1991,  Momot  1992). 

Watersheds:  2HG  &  2KD  -  Scugog  River  &  Madawaska  River 

Orconectes  virilis  was  the  only  crayfish  species  caught  in  Scugog  Lake  (Table  8). 
Similarly,  C.  bartoni  was  the  only  crayfish  found  in  Timberwolf  Lake.  Both  of  these 
species  were  caught  in  Louisa  Lake.  No  other  crayfish  species  were  found  in  these  2 
watersheds,  although  more  species  will  probably  be  collected  when  other  lakes  are 
surveyed.  For  example,  two  of  the  14  temporal  reference  lakes  (see  David  et  al.  1994) 
are  in  the  Upper  Madawaska  watershed.  One  of  these  lakes,  Cradle  Lake,  only 
contains  a  population  of  C.  bartoni.  By  contrast  a  second  lake,  Delano  Lake,  supports 
populations  of  O.  virilis.  O.  propinquus  and  C.  bartoni.  Similarly,  Somers  and  Green 
(1993)  caught  O.  obscurus  in  Rock  Lake,  which  is  downstream  of  Louisa  Lake  in 
Algonquin  Provincial  Park.  Until  more  lakes  have  been  surveyed  in  these  two 
watersheds  it  is  not  appropriate  to  compare  the  relative  frequencies  of  occurrence  of 
the  various  crayfish  species  with  data  from  the  other  watersheds. 

Environmental  Correlates  with  Crayfish  Relative  Abundance 

Simple  correlations  between  total  crayfish  CUE  and  the  environmental  variables 
indicated  no  strong  relationships  (see  Table  9).  Scatterplots  revealed  several  unusually 
high  catches  generally  attributed  to  a  single  species  (e.g.,  the  catch  of  more  than  12  C. 
bartoni  per  trap  in  Clayton  Lake).  Once  these  outliers  were  removed,  linear  or  even 
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curvilinear  (i.e.,  rank-order)  patterns  were  still  weak  (e.g.,  see  Figure  3).  Further 
examination  of  the  scatterplots  revealed  that  each  of  the  3  common  species  (i.e.,  C. 
barton!,  O.  propinquus  and  O.  virilis)  displayed  unique  patterns  associated  with  different 
environmental  variables. 

Both  Pearson  and  Spearman  rank  correlations  between  the  trap  catches  of  the 
individual  crayfish  species  and  the  environmental  variables  were  also  weak  (Table  9). 
For  example,  only  the  rank  correlations  between  the  CUE  of  C.  barton i  and  manganese 
concentration  exceeded  0.30.  Somewhat  higher  values  were  observed  for  the  CUE  of 
C.  robustus.  although  these  correlations  were  only  based  on  data  from  12  lakes  (Table 
8).  By  comparison,  Pearson  correlations  for  the  catch  of  O.  propinquus  were  all  less 
than  0.30,  although  several  rank  correlations  exceeded  0.45.  Similarly  there  were  no 
strong  environmental  correlations  with  the  catches  of  O.  virilis.  These  results  suggest 
that  linear  and  monotonic  models  of  the  relationships  between  crayfish  CUE  and  the 
environmental  variables  are  unrealistic. 

A  subset  of  the  relationships  between  crayfish  relative  abundance  and  the 
environmental  variables  are  illustrated  with  plots  of  the  CUE  of  the  3  common  species 
and  elevation,  pH  and  aluminum  concentration  (Figure  3).  Elevation  may  influence  the 
distribution  of  crayfish  within  watersheds  because  high-elevation  lakes  are  smaller  on 
average,  with  less  glacial  till  in  their  watersheds  (Jefferies  and  Snyder  1983,  Girard  and 
Reid  1990).  As  a  result,  these  lakes  tend  to  be  more  susceptible  to  lake  acidification 
because  of  reduced  buffering  capacity  and  smaller  lake  volumes  (Dillon  et  al.  1984, 
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1987,  France  and  Collins  1993).  Moreover,  lakes  at  higher  elevations  nriay  be  more 
inaccessible  for  crayfish  that  dispersed  from  glacial  refugia  (Crocker  and  Barr  1968, 
Guiasuetal.  1996). 

David  et  al.  (1994)  found  evidence  suggesting  that  declines  in  crayfish  CUE  coincided 
with  increases  in  aluminum  concentrations  in  Dorset-area  lakes  that  were  monitored  for 
five  years.  Laboratory  toxicology  studies  (e.g.,  Malley  1980,  Berrill  et  al.  1985,  Malley 
and  Chang  1985)  have  also  found  that  low  pH  and  high  aluminum  are  toxic  to  several 
orconectid  crayfish  species.  Thus,  pH  and  aluminum  may  also  be  important  factors 
affecting  the  abundance  and  distribution  of  crayfish  in  south-central  Ontario  lakes  (e.g., 
see  France  1993,  France  and  Collins  1993). 

From  these  simple  scatterplots,  we  found  that  both  C.  bartoni  and  O.  virilis  were  most 
abundant  in  high-elevation  lakes  (Figure  3).  This  observation  suggests  that  factors 
other  than  lake  acidification  may  reduce  the  abundances  of  these  two  species  in  low- 
elevation  lakes.  For  example,  crayfish  catches  in  low-elevation  lakes  may  be  affected 
by  greater  numbers  of  predatory  fish  species  (especially  bass;  Mather  and  Stein  1993, 
Hinch  et  al.  1994),  or  perhaps  by  competition  with  other  crayfish  species  (Olsen  et  al. 
1 991 ,  Hill  and  Lodge  1 994,  Guiasu  et  al.  1 996).  By  contrast,  O.  propinquus  was 
abundant  in  lakes  across  the  range  of  elevations.  These  differences  in  CUE  associated 
with  elevation  emphasize  that  the  individual  species  illustrate  different  responses  to 
geographic  features  such  as  elevation,  and  hence,  that  pooling  catch  data  for  all 
species  will  lose  resolution  by  combining  different  (and  potentially  conflicting)  patterns. 

43 


The  plots  of  CUE  with  pH  provide  a  scatter  of  points  that  is  similar  to  the  plots  with 
elevation  (Figure  3).  C.  bartoni  was  collected  from  lakes  spanning  a  broad  range  of  pH 
(I.e.,  from  5.0  to  >7.0)  with  an  optimum  or  peak  abundance  in  lakes  with  a  pH  of 
approximately  6.5.  Orconectes  propinquus  was  also  found  in  lakes  spanning  a  broad 
pH  range;  however,  O.  propinquus  was  caught  in  only  two  lakes  with  a  pH  below  5.8 
(i.e..  Round  and  Black  lakes  in  watershed  2EA  with  pH  values  of  5.0  and  4.9, 
respectively;  Table  2).  Unlike  C.  bartoni.  the  CUE  of  O.  propinquus  failed  to  illustrate  a 
clear  optimum  or  peak  abundance,  although  the  catches  approached  a  maximum  in 
lakes  with  a  pH  near  6.7.  By  contrast,  no  O.  virilis  were  caught  in  lakes  with  a  pH  below 
5.7.  This  pattern  is  consistent  with  observations  on  the  demise  of  a  population  of  O. 
virilis  in  a  lake  that  was  experimentally  acidified  in  northwestern  Ontario  (France  1987, 
Davies  1989). 

Each  species  illustrates  somewhat  similar  patterns  with  respect  to  aluminum 
concentration  (Figure  3).  For  example,  the  largest  catches  of  C.  bartoni  were  from 
lakes  with  aluminum  concentrations  below  60  ^g-L'^  (with  one  exception.  Round  Lake  in 
watershed  2EA  with  105  A^g-L"^).  Catches  of  O.  propinquus  were  also  higher  in  lakes 
with  aluminum  levels  below  approximately  60  lu.g-L'^  (again  with  one  exception.  Round 
Lake).  Similarly,  O.  virilis  abundances  approached  zero  in  lakes  where  the  aluminum 
concentrations  exceeded  50-75  lug-L'^  (with  two  marked  exceptions,  Fogal  and  Fawn  2 
lakes  in  watershed  2EB). 

These  trends  in  CUE  with  respect  to  3  environmental  variables  illustrate  that  each 
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crayfish  species  responds  to  environmental  factors  in  somewhat  different  ways.  As  a 
result,  a  composite  measure  such  as  the  total  crayfish  catch  will  obscure  the  individual 
relationships  between  crayfish  relative  abundance  and  the  various  environmental 
factors.  Moreover,  the  graphs  reveal  that  the  relationships  between  crayfish  CUE  and 
the  environmental  variables  are  not  simple  linear  or  curvilinear  functions.  As  a  result,  a 
simple  correlation  matrix,  or  even  a  rank-correlation  matrix,  will  not  adequately 
summarize  the  environmental  patterns  in  crayfish  relative  abundance. 

Because  of  these  problems,  we  also  used  a  multiple  regression  with  linear,  squared 
and  cubed  terms  (i.e.,  a  polynomial  regression)  to  examine  curvilinear  relationships 
between  crayfish  CUE  and  the  22  environmental  parameters  (Table  9).  The  square 
root  of  the  resultant  multiple  correlation  is  an  index  of  the  predictive  power  of  the 
polynomial  regression  that  can  be  compared  with  the  Pearson  and  Spearman 
correlations.  On  average,  the  multiple  correlations  were  noticeably  higher  than  the 
Pearson  or  Spearman  rank  correlations.  However,  none  of  the  multiple  correlations 
using  total  crayfish  CUE  was  larger  than  0.225.  This  finding  is  consistent  with  our 
proposal  that  total  crayfish  CUE  obscures  the  separate  patterns  associated  with  each 
species. 

The  tabulated  correlations  from  the  polynomial  regression  reveal  one  additional 
problem  (Table  9).  Almost  all  of  the  multiple  correlations  for  C.  robustus  exceed  0.30. 
This  undoubtedly  occurs  because  of  the  relatively  small  number  of  lakes  that  were  used 
in  these  regressions  (i.e.,  C.  robustus  was  only  captured  in  12  lakes).  Consequently, 

45 


we  advise  caution  when  examining  the  results  for  C.  robustus  because  the  number  of 
parameters  in  the  regression  model  is  large  relative  to  the  number  of  observations.  By 
contrast,  few  of  the  multiple  correlations  for  the  other  taxa  are  very  large.  For  example, 
the  3  largest  values  for  C.  bartoni  involve  magnesium,  conductivity  and  maximum 
depth.  The  3  largest  multiple  correlations  with  the  CUE  of  O.  propinquus  involve 
maximum  depth,  NH4  and  aluminum  concentration.  Similarly  maximum  depth, 
elevation,  plus  Na,  TKN  and  DOC  exhibited  the  largest  correlations  with  the  CUE  of  O. 
virilis.  Although  maximum  depth  was  an  important  predictor  for  each  species,  no  other 
variable  was  highly  correlated  with  the  CUE  of  each  crayfish  species.  These  results 
suggest  that  a  simple  regression  approach  to  model  crayfish  CUE  is  not  very  effective, 
and  that  the  polynomial  regression  only  provides  a  modest  improvement. 

In  addition  to  the  regression  analysis  using  the  individual  species,  we  examined  the 
relationships  between  multivariate  trends  in  crayfish  CUE  and  the  environmental 
variables.  Here  we  capitalized  on  the  concepts  behind  canonical  correspondence 
analysis  (CCA)  which  utilizes  a  simple  correspondence  analysis  (CA)  ordination 
followed  by  multiple  regression  to  determine  the  environmental  variables  that  best 
predicts  the  CA  scores  (ter  Braak  1988,  1990,  Palmer  1993).  However,  here  we  used 
separate  polynomial  regressions  between  the  CA  axes  and  each  of  the  22 
environmental  variables  because  of  relatively  high  multicollinearity  among  the 
morphology  and  water  chemistry  variables.  Analyses  of  simulated  and  field  data  have 
illustrated  that  CA  is  a  relatively  robust  ordination  technique  that  is  suitable  for  species 
abundance  data  that  span  environmental  gradients  (Digby  and  Kempton  1987,  ter 
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Braak  and  Prentice  1988).  Given  the  weak  correlations  from  the  raw  CUE  data  (Table 
9),  the  CA  approach  has  appeal  because  it  provides  a  multivariate  summary  of  the 
major  (nonlinear)  trends  in  the  crayfish  CUE  data. 

The  CA  ordination  of  the  log-transformed  crayfish  CUE  summarized  71.5%  of  the 
overall  variation  in  the  first  two  axes  (i.e.,  CA1  and  CA2).  The  third  axis  incorporated 
another  20.4%  (i.e.,  CAS),  resulting  in  92%  of  the  variation  in  the  first  3  axes.  Only  the 
71  lakes  where  crayfish  were  captured  were  used  in  the  ordination.  The  relative 
positions  of  the  lakes  are  indicated  in  the  plot  of  the  first  two  axes  (Figure  4).  Lakes 
that  are  close  together  in  the  ordination  have  similar  species  compositions  and  similar 
relative  abundances.  For  example.  Little  Margaret  (LM)  and  Clear  (CA)  lakes  are 
neighbours  in  the  upper  left  quadrant  of  the  plot.  Cambarus  bartoni  and  O.  virilis  were 
caught  in  both  lakes.  The  catches  of  C.  bartoni  were  0.09  and  0.06  crayfish  per  trap, 
respectively;  whereas  the  catches  of  O.  virilis  were  3.28  and  2.83,  respectively  (see 
Table  8). 

The  species  vectors  were  also  plotted  to  assist  with  the  interpretation  of  the  lake  scores 
(Figure  4).  The  3  common  taxa,  C.  bartoni.  O.  propinquus  and  O.  virilis,  were  the  most 
important  taxa  with  respect  to  the  first  two  CA  axes.  The  species  vectors  indicate  that 
the  abundance  of  C.  bartoni  increases  from  the  origin  towards  the  upper  right  of  the 
figure.  Orconectes  virilis  increases  towards  the  upper  left,  whereas  O.  propinquus 
increases  in  abundance  towards  the  lower  left.  By  reference  to  the  plot  of  the  lake 
scores,  lakes  lying  along  any  one  of  these  3  vectors  have  catches  predominated  by  that 
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species.  Lakes  lying  between  any  two  of  these  vectors  have  catches  comprised  of  both 
species.  For  example,  the  vector  represented  by  5  lakes  that  extends  from  the  central 
cluster  of  lakes  beginning  with  GS  and  ending  at  CY  spans  a  gradient  of  increasing 
catches  of  C.  bartoni  from  Grass  (0.52  crayfish  per  trap  night;  Table  8)  to  Clayton  Lake 
(12.78  per  trap  night).  Similarly,  Grandview  (GV)  and  Bigwind  (BW)  lakes  had  large 
catches  of  O.  propinquus  (i.e.,  3.07  and  8.68  per  trap  night,  respectively)  with  small 
catches  of  C.  bartoni  in  Grandview  (0.06)  and  O.  virilis  in  Bigwind  (1 .26).  These  smaller 
catches  of  a  second  species  pull  the  lakes  away  from  the  O.  propinquus  vector  towards 
the  C.  bartoni  or  O.  virilis  vectors.  In  a  similar  context,  Ridout  (RT)  and  Little  Margaret 
(LM)  lakes  had  the  largest  catches  of  O.  virilis  (i.e.,  4.35  and  3.28,  respectively).  These 
two  lakes  occupy  extreme  positions  in  the  upper  left  quadrant. 

The  first  3  CA  axes  were  correlated  with  the  22  environmental  variables  to  search  for 
linear,  monotonic  and  polynomial  trends  (Table  10).  Only  3  Pearson  correlations  were 
greater  than  0.30.  These  correlations  involved  CA3  and  SO4  (r=-0.516),  Mn  (r=-0.473) 
and  pH  (r=0.325).  Similarly,  none  of  the  Spearman  rank  correlations  with  the  first  or 
second  CA  axes  were  larger  than  0.30.  The  largest  rank  correlations  with  the  third  CA 
axis  Included  alkalinity,  TKN,  DIC,  pH  and  elevation  (r3=0.399,  0.335,  0.332,  0.327,  and 
-0.327,  respectively).  These  trends  suggest  that  lakes  separated  along  CA3  span  a 
range  from  high  alkalinity,  high  pH  and  low  elevation  at  the  positive  end  of  the  axis,  to 
low  alkalinity,  low  pH,  high  elevation  and  high  sulphate  concentrations  at  the  negative 
end  of  the  axis. 
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Most  of  the  multiple  correlations  from  the  polynomial  regressions  were  also  less  than 
0.30  (Table  10).  Only  the  correlation  between  CA1  and  NH4  exceeded  0.30,  although 
the  correlation  with  Al  approached  this  value  (i.e.,  0.294).  Two  environmental  variables 
were  correlated  with  the  second  CA  axis  (i.e.,  DOC  and  NH4).  Although  the  first  two  CA 
axes  summarize  patterns  that  can  be  attributed  to  the  relative  abundances  of  C. 
bartoni.  O.  propinquus  and  O.  virilis.  variation  represented  by  these  two  axes  was  only 
weakly  correlated  with  the  environmental  variables.  Perhaps  some  other  factor,  such 
as  the  number  and  types  of  predatory  fish  species,  contributes  to  the  observed  pattern 
in  the  lake  scores  on  the  first  two  CA  axes. 

By  contrast,  more  of  the  environmental  variables  were  correlated  with  CA3.  For 
example,  the  polynomial  regressions  for  6  of  the  22  variables  produced  correlations  that 
exceeded  0.30.  These  correlations  included  SO4  (0.627),  Mn  (0.611),  DOC  (0.401), 
alkalinity  (0.383),  DIC  (0.365)  and  pH  (0.337).  In  a  biological  context,  lakes  with 
negative  scores  on  CA3  support  populations  of  C.  robustus  (and  to  a  lesser  extent,  O. 
immunis  and  O.  obscurus),  whereas  lakes  with  positive  scores  support  populations  of 
O.  propinquus  or  C.  bartoni.  These  results  reflect  the  analysis  of  the  individual  species 
CUEs  (e.g.,  Figure  3),  since  lakes  containing  O.  propinquus  were  characterized  by  high 
pH,  alkalinity,  DIC  and  TKN,  and  low  elevation,  plus  low  sulphate  and  manganese 
concentrations.  These  patterns  can  be  confirmed  by  comparing  the  water  chemistry  of 
Koshlong  (KG)  and  Grandview  (GV)  lakes  which  lie  at  opposite  ends  of  CA3  (Table  2). 

This  multivariate  analysis  was  repeated  using  crayfish  presence-absence  data  for  the 
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71  lakes  where  crayfish  were  caught.  The  first  two  CA  axes  summarized  almost  62%  of 
the  variation  in  the  resultant  ordination  (i.e.,  39.0  and  22.8%  in  CM  and  CA2, 
respectively;  Figure  5).  The  third  CA  axis  incorporated  another  17.3%  of  the  variation 
(i.e.,  CA3).  Despite  using  71  lakes  in  the  CA,  only  20  unique  sets  of  lake  scores  were 
obtained  on  the  first  two  CA  axes.  The  resultant  plot  revealed  that  each  of  these 
coordinates  represents  a  particular  combination  of  crayfish  species.  For  example,  lakes 
with  C.  bartoni  lie  to  the  right  of  zero  on  CA1 .  Lakes  where  O.  propinquus  was  caught 
lie  above  a  line  at  an  angle  of  about  15  degrees  to  CA1 ,  whereas  lakes  with  O.  virilis  fall 
below  a  line  lying  at  an  angle  of  about  165  degrees  to  CA1 .  An  inset  of  the  species 
scores  on  the  first  two  CA  axes  supports  this  interpretation  (Figure  5). 

Pearson  correlations  between  the  CA  scores  from  the  crayfish  presence-absence 
matrix  and  the  22  environmental  variables  revealed  trends  that  differed  from  the  CA  of 
the  logarithmically  transformed  abundance  data  (Table  11).  Only  4  Pearson 
correlations  were  larger  than  0.30,  and  all  of  these  were  associated  with  CA1  (i.e., 
elevation,  pH,  Fe  and  DOC;  0.490,  -0.353,  -0.325  and  -0.322,  respectively).  In  addition, 
correlations  between  CA1  and  maximum  depth  (r=0.299)  and  total  phosphorus 
concentration  (r=-0.297)  approached  the  0.30  threshold.  By  contrast,  all  of  the  Pearson 
correlations  with  CA2  or  CA3  were  less  than  0.30. 

Many  of  the  Spearman  rank  correlations  were  larger  than  the  corresponding  Pearson 
correlations,  again  suggesting  that  the  relationships  between  the  CA  axes  and  the 
environmental  variables  were  curvilinear  (Table  1 1 ).  For  example,  1 3  of  the  22  rank 
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correlations  with  CA1  exceeded  0.30.  These  correlations  included  alkalinity  (-0.503), 
Die  (-0.502),  elevation  (0.468),  Na  (-0.406),  pH  (-0.401),  TP  (-0.397),  maximum  depth 
(0.389),  DOC  (-0.377),  Cond25  (-0.362),  TKN  (-0.357),  Mg  (-0.346),  Ca  (-0.336)  and  CI 
(-0.306).  Like  the  Pearson  correlations,  none  of  the  rank  correlations  with  CA2  or  CA3 
exceeded  0.30. 

Results  from  the  polynomial  regressions  confirmed  the  curvilinear  nature  of  the 
relationships  between  the  presence-absence  CA  scores  and  the  environmental 
variables  (Table  11).  Half  of  the  correlations  between  CA1  and  the  22  environmental 
variables  were  larger  than  0.30.  These  correlations  included  elevation  (0.490),  DIC 
(0.454),  alkalinity  (0.439),  maximum  depth  (0.393),  TP  (0.384),  DOC  (0.370),  pH 
(0.356),  Fe  (0.343),  SO^  (0.338),  TKN  (0.324)  and  Mg  (0.302).  Only  3  correlations  with 
CA2  (i.e.,  K  [0.481],  Mg  [0.318]  and  SO4  [0.300])  and  none  of  the  multiple  correlations 
with  CA3  exceeded  0.30. 

The  multiple  correlations  between  CA1  and  the  22  environmental  variables  suggest  that 
the  scatter  of  lakes  along  CA1  reflects  aspects  of  the  water  chemistry  and  lake 
morphology.  Lakes  lying  at  the  positive  end  of  the  first  axis  are  at  higher  elevation,  with 
lower  pH,  lower  alkalinity  and  lower  conductivity.  Lakes  at  the  negative  end  of  the  axis 
will  be  lower  in  elevation  with  relatively  shallower  maximum  depths  and  higher 
concentrations  of  Ca,  Mg  and  nutrients. 

Because  only  one  axis  from  each  of  the  two  ordinations  exhibited  nontrivial  correlations 
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with  the  set  of  environmental  variables,  we  plotted  the  lake  scores  for  these  two  axes  to 
establish  whether  the  underlying  patterns  were  similar  (Figure  6).  The  resultant 
Pearson  correlation  was  0.194,  suggesting  that  the  two  axes  summarized  different 
trends.  For  example,  lakes  with  positive  scores  on  the  presence-absence  data  CA1  are 
lakes  with  populations  of  C.  bartoni.  Lakes  with  positive  scores  on  the  log-abundance 
CAS  tend  to  have  higher  abundances  of  C.  bartoni  and,  or  O.  propinquus.  As  a  result, 
this  CAS  separates  lakes  with  high  abundances  of  these  two  species  from  lakes  with 
other  species  combinations.  This  interpretation  is  supported  by  the  fact  that  Clayton 
(CY),  Grandview  (GV),  Red  Chalk  East  (RE)  and  Bigwind  (BW)  had  the  largest  catches 
of  these  two  species. 

Lakes  lying  at  the  negative  end  of  CAS  had  relatively  high  abundances  of  C.  robustus 
(Figure  6).  This  trend  of  particular  species,  or  pairs  of  species  contributing  to  the 
separation  of  subsets  of  lakes  suggested  that  perhaps  the  presence  or  absence  of 
species  (e.g.,  see  Guiasu  et  al.  1996),  rather  than  their  relative  abundances  was  more 
readily  predicted  from  the  environmental  data.  This  proposal  is  further  supported  by  the 
fact  that  many  of  the  environmental  variables  were  correlated  with  the  first  axis  of  the 
presence-absence  CA  (i.e.,  the  most  important  axis),  as  compared  to  the  third  CA  axis 
from  the  log-abundance  data.  We  examined  this  hypothesis  using  one-way  analysis  of 
variance  (ANOVA)  to  compare  the  mean  values  for  the  various  environmental  variables 
for  lakes  with  versus  lakes  without  various  crayfish  species. 

The  first  series  of  ANOVAs  compared  the  means  of  the  environmental  variables  for 
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lakes  where  crayfish  were  captured  (N=71)  to  the  means  for  the  29  lakes  where  no 
crayfish  were  caught  (Table  12).  Only  two  of  the  22  ANOVAs  revealed  significant 
differences  (P<0.05;  i.e.,  the  mean  concentrations  of  NO3  and  SOJ.  In  addition, 
differences  in  elevation  approached  significance  (i.e.,  0.10>P>0.05).  Crayfish  were 
absent  from  the  higher-elevation  lakes,  or  lakes  with  high  concentrations  of  NO3  or  SO4. 
By  contrast,  the  multivariate  ANOVA  (MANOVA)  incorporating  all  22  variables  was  not 
significant  (F=1.166  with  22  and  77  df,  P=0.303),  suggesting  that  the  simple  presence 
or  absence  of  crayfish  could  not  be  predicted  using  the  22  variables  simultaneously. 

This  analysis  was  extended  to  contrast  lakes  on  the  basis  of  the  number  of  crayfish 
species  captured.  Again,  one-way  ANOVAs  were  used  to  compare  the  mean  values  for 
each  of  the  22  environmental  variables  for  lakes  with  0  (N=29),  1  (N=23),  2  (N=34)  and 
3-or-more  species  (N=14;  Table  12).  The  ANOVAs  revealed  that  differences  in 
elevation  and  NO3  concentration  approached  significance  (0.10>P>0.05),  but  none  of 
the  other  variables  was  significantly  different  (P>0.10).  As  above,  these  one-way 
ANOVAs  re-iterated  that  lakes  without  crayfish  tended  to  be  at  higher  elevation  with 
higher  concentrations  of  NO3.  The  absence  of  significant  differences  in  the  univariate 
comparisons  was  confirmed  by  the  corresponding  MANOVA  (F=1.029  with  66  and  225 
df,  P=0.429).  Consequently,  none  of  the  22  environmental  variables,  alone  or  in 
combination,  could  adequately  predict  the  number  of  crayfish  species  caught  in  a  given 
lake.  This  observation  is  consistent  with  earlier  findings  that  total  abundances  were 
only  weakly  correlated  with  the  environmental  variables  because  each  species 
displayed  somewhat  different  patterns  across  the  various  gradients. 
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One-way  ANOVAs  comparing  mean  values  for  the  22  environmental  variables  based 
on  the  presence  or  absence  of  the  3  common  crayfish  species  indicated  that  different 
environmental  variables  were  associated  with  the  presence  or  absence  of  each  species 
(Tables  13  and  14).  For  example,  the  mean  values  for  lakes  supporting  populations  of 
O.  propinquus  were  significantly  different  from  the  means  for  lakes  without  O. 
propinquus  with  respect  to  elevation  and  pH  (P<0.05).  In  addition,  differences  in  lake 
maximum  depth,  DOC,  SO4  and  TP  concentrations  approached  significance 
(0.10>P>0.05).  The  associated  MANOVA  was  significant  (F=1. 717  with  22  and  77  df, 
P<0.044),  indicating  that  lakes  with  O.  propinquus  could  be  distinguished  from  lakes 
where  O.  propinquus  was  absent.  In  general,  lakes  supporting  populations  of  O. 
propinquus  were  at  lower  elevation,  with  higher  pH,  shallower  maximum  depths,  and 
higher  concentrations  of  DOC  and  TP,  but  lower  concentrations  of  SO4. 

Somewhat  different  variables  were  associated  with  the  presence  or  absence  of  O.  virilis 
(Table  13).  For  example,  differences  in  lake  elevation  were  relatively  unimportant 
(P=0.198).  However,  significant  differences  in  average  pH,  TP  and  Mn  concentrations 
were  observed  (P<0.05).  In  addition,  differences  in  the  mean  concentration  of  Al,  plus 
conductivity  and  alkalinity  approached  significance  (i.e.,  0.10>P>0.05).  In  general, 
lakes  supporting  O.  virilis  populations  were  characterized  by  higher  pH,  higher  TP  and 
lower  Mn  concentrations.  Alkalinity  and  DIC  concentrations  were  also  lower,  and  Al 
concentrations  were  higher  in  lakes  where  O.  virilis  were  absent.  Although  several  of 
the  univariate  comparisons  were  significant,  the  MANOVA  only  approached  significance 
(F=1.656  with  22  and  77  df,  P=0.055).  Consequently,  the  environmental  gradient 
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associated  with  the  presence  or  absence  of  O.  virilis  was  somewhat  less  pronounced 
than  that  observed  for  O.  propinquus. 

Mean  values  for  4  of  the  22  environmental  variables  were  significantly  different  between 
lakes  with  populations  of  C.  bartoni  versus  lakes  where  C.  bartoni  was  absent  (Table 
14).  Cambarus  bartoni  was  caught  in  the  deeper  lakes  at  higher  elevations  (P<0.05). 
Chemically,  these  lakes  were  characterized  by  lower  Fe  and  DOC  concentrations 
(P<0.05),  plus  lower  TP,  TKN,  Cond25,  pH  and  DIC  concentrations  (0.10>P>0.05). 
The  MANOVA  also  revealed  significant  differences  between  the  lakes  with  and  without 
C.  bartoni  (F=1.832  with  22  and  77  df,  P=0.028). 

A  simple  comparison  of  the  multivariate  F  values  for  these  3  species  indicates  that  the 
occurrence  of  C.  bartoni  could  be  predicted  with  better  precision  than  O.  propinquus  or 
O.  virilis.  Cambarus  bartoni  was  found  in  deep,  high-elevation  lakes  that  were 
somewhat  acidic  with  low  DOC  concentrations.  This  finding  Is  consistent  with  similar 
studies  of  lakes  in  south-central  Ontario  (e.g.,  Reid  and  David  1990,  David  et  al.  1994). 
Laboratory  toxicological  studies  by  Berrill  et  al.  (1985)  have  suggested  that  C.  bartoni  is 
relatively  tolerant  of  low-pH  conditions  when  compared  to  several  orconectids.  Since 
the  smaller,  high-elevation  lakes  in  south-central  Ontario  are  most  vulnerable  to  the 
historical  impacts  of  acid  precipitation  (Dillon  et  al.  1984,  1987,  France  and  Collins 
1993),  it  is  not  surprising  that  an  intolerant  species  like  C.  bartoni  is  most  common  in 
these  lakes.  Similarly,  the  acid-sensitive  species  O.  propinquus  was  infrequently  found 
in  high-elevation  lakes. 
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Unfortunately,  the  distribution  of  O.  virilis  was  difficult  to  resolve.  Lake  elevation  was 
not  important  (P>0.05,  Table  13),  but  pH,  alkalinity,  Al  and  TP  concentrations  were. 
These  water-chemistry  parameters  are  indicators  of  lake  acidification  (Dillon  et  al. 
1984).  Malley  (1980),  Malley  and  Chang  (1985),  France  (1987),  Davies  (1989)  and 
France  and  Collins  (1993)  all  document  the  sensitivity  of  O.  virilis  to  changes  in  water 
chemistry  associated  with  lake  acidification.  In  all  likelihood,  the  distribution  of  O.  virilis 
reflects  a  response  to  the  historical  acidification  of  some  of  the  survey  lakes  (e.g., 
France  and  Collins  1993).  However,  this  pattern  is  undoubtedly  confounded  by 
recolonization  (e.g.,  Gunn  and  Keller  1990)  and  other  biogeographic  factors  including 
the  distribution  of  predatory  fish  (e.g.,  Emery  1975,  Capelli  and  Magnuson  1983,  Collins 
et  al.  1983,  DiDonato  and  Lodge  1993).  More-careful  scrutiny  of  sets  of  lakes  within 
watersheds  may  assist  in  resolving  the  confusing  effects  of  elevation,  water  chemistry 
and  fish  predation. 

Recommendations 

One  objective  of  this  survey  was  to  develop  predictive  models  of  crayfish  distribution 
patterns  in  south-central  Ontario  watersheds.  The  data  for  4  of  the  7  watersheds  are 
fairly  comprehensive,  but  more  sampling  needs  to  be  done  in  watersheds  2CF,  2HG 
and  2KD.  Much  like  the  temporal  crayfish  study  (David  et  al.  1994),  a  standard 
sampling  protocol  must  be  followed  to  ensure  consistency  for  both  sampling  and  data 
processing.  Consequently,  more  lakes  in  the  2CF,  2HG  and  2KD  watersheds  should 
be  surveyed  using  the  sampling  protocol  described  in  this  report  (also  see  Reid  and 
David  1990,  David  etal.  1994). 
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Many  factors  must  be  considered  when  comparing  crayfish  catches  from  a  group  of 
lakes.  Finding  crayfish  Is  contingent  upon  parameters  such  as  life  history  stages  (e.g., 
date  of  sampling:  Somers  and  Green  1993),  water  chemistry  (e.g.,  Schindler  1987, 
France  1993),  habitat  (Jones  and  Momot  1981),  the  presence  of  predatory  fish  (Collins 
et  al.  1983,  Somers  and  Green  1993),  dispersal  (Crocker  and  Barr  1968,  Capelli  and 
Magnuson  1983,  Guiasu  et  al.  1996),  and  the  presence  of  other  crayfish  species 
including  introduced  species  (e.g.,  Berrill  1978,  Momot  1992,  DiDonato  and  Lodge 
1993,  Kershner  and  Lodge  1995,  Guiasu  et  al.  1996).  Herein  we  have  found  that 
broad-scale  surveys  that  span  several  watersheds  will  not  necessarily  resolve  complex 
biogeographic  patterns  (Hinch  et  al.  [1994]  report  similar  problems).  In  all  likelihood, 
better  resolution  would  be  obtained  by  focusing  on  a  larger  number  of  lakes  within  one- 
or-two  watersheds.  This  type  of  approach  should  span  an  elevation  gradient  with  a 
minimum  of  confounding  factors. 

This  survey  encountered  two  crayfish  species  that  have  been  introduced  into  lakes  in 
south-central  Ontario  (i.e.,  O.  obscurus  and  O.  rusticus  were  found  in  2  of  the  7 
watersheds).  Although  the  research  on  O.  obscurus  is  limited  (e.g.,  Crocker  and  Barr 
1968,  Berrill  1978,  Corey  1988),  considerably  more  information  is  available  for  O. 
rusticus.  For  example,  this  latter  species  is  highly  competitive  and  may  displace  native 
species  (e.g.,  O.  propinquus  and  O.  virilis:  Lodge  et  al.  1986,  Olsen  et  al.  1991,  Hill  and 
Lodge  1994).  Recently,  O.  rusticus  has  been  reported  in  the  Lake  Superior  basin  in 
northwestern  Ontario  (Momot  1992).  In  addition  to  displacing  native  species,  Momot 
(1992)  suggests  that  O.  rusticus  may  also  alter  the  benthic  invertebrate  community.   In 
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order  to  understand  these  types  of  impacts,  complementary  data  on  crayfish  catches 
and  benthic  invertebrate  collections  (e.g.,  David  et  al.  1994,  Reid  et  al.  1994  a,  b) 
should  be  compared  to  search  for  potential  relationships  that  may  change  as  crayfish 
species  assemblages  change. 

Since  O.  rusticus  is  so  competitive  and  has  been  progressively  expanding  its  range,  it 
would  be  advisable  to  monitor  lakes  wherever  this  species  is  found.  This  may  require  a 
more  frequent  sampling  schedule  than  those  lakes  where  introduced  species  have  not 
been  found.  For  example,  the  temporal  crayfish  study  (David  et  al.  1994)  revealed 
significant  declines  in  crayfish  abundances  over  a  5-year  period  in  several  south-central 
Ontario  lakes.  Consequently,  we  propose  that  those  lakes  where  introduced  crayfish 
species  have  been  found  should  be  sampled  at  least  once  every  5  years. 

The  preliminary  QA/QC  results  from  sampling  7  lakes  on  two  separate  dates  indicate  a 
need  for  routine  QA/QC  checks  (also  see  Reid  et  al.  1995).  At  a  minimum,  several 
lakes  should  be  resampled  during  a  given  field  season  to  provide  better  estimates  of 
the  magnitude  of  within-season  sampling  variablity.  Ultimately  this  information,  when 
combined  with  estimates  of  within-lake  temporal  variation,  should  be  used  to  develop 
estimates  of  the  detection  limits  for  the  standard  crayfish-sampling  protocol.  These 
detection  limits  are  critical  in  the  assessment  of  presence-absence  data  and  the 
evaluation  of  temporal  changes  in  crayfish  abundances. 

Lastly,  we  propose  that  year-to-year  changes  in  the  sex  ratios  of  crayfish  catches  may 
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reflect  the  influence  of  global  warming  on  seasonal  patterns  in  the  annual  warming  and 
cooling  cycles  in  small,  inland  lakes.  This  relationship  should  be  examined  as 
additional  data  are  added  to  the  temporal  data  series. 
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Table  3:  Fish  species  names  and  OMNR  code  numbers 


OMNR  Code 

Scientific  Name 

Common  Name 

76 

Oncorhvnchus  mvkiss 

Rainbow  Trout 

80 

Salvelinus  fontinalis 

Brook  Trout 

81 

Salvelinus  namavcush 

Lake  Trout 

82 

Salvelinus  namavcush  x  fontinalis  Solake 

91 

Coreqonus  cluDeafonmis 

Lake  Whitefish 

93 

Coreqonus  artedii 

Cisco 

121 

Osmerus  mordax 

Rainbovi^  Smelt 

131 

Esox  lucius 

Northern  Pike 

132 

Esox  masQuinongv 

Muskellunge 

141 

Umbra  limi 

Central  Mudminnow 

162 

Catostomus  catostomus 

Longnose  Sucker 

163 

Catostomus  commersoni 

White  Sucker 

182 

Phoxinus  eos 

Northern  Redbelly  Dace 

185 

Couesius  Diumbeus 

Lake  Chub 

186 

CvDrinus  carpio 

Common  Carp 

189 

Hvboqnathus  hankinsoni 

Brassy  Minnow 

194 

Notemiaonus  crvsoleucas 

Golden  Shiner 

196 

NotroDis  athennoides 

Emerald  Shiner 

198 

NotroDis  comutus 

Common  Shiner 

199 

NotroDis  heterodon 

Blackchin  Shiner 

200 

Notropis  heteroleois 

Blacknose  Shiner 

201 

NotroDis  hudsonius 

Spottail  Shiner 

206 

Notropis  volucellus 

Mimic  Shiner 

208 

Pimephales  notatus 

Bluntnose  Minnow 

209 

Pimeohales  promelas 

Fathead  Minnow 

210 

Rhinichthvs  atratulus 

Blacknose  Dace 

212 

Semotilus  atromaculatus 

Creek  Chub 

214 

Semotilus  maraarita 

Pearl  Dace 

232 

Ictalurus  natalis 

Yellow  Bullhead 

233 

Ictalurus  nebulosus 

Brown  Bullhead 

261 

Fundulus  diaphanus 

Banded  Killifish 

271 

Lota  lota 

Burbot  (ling) 

281 

Culaea  inconstans 

Brook  Stickleback 

291 

Percopsis  omiscomavcus 

Trout-perch 

311 

Ambloplites  rupestris 

Rock  Bass 

313 

Lepomis  qibbosus 

Pumpkinseed 

314 

Lepomis  macrochirus 

Bluegill 

316 

Microoterus  dolomieui 

Smallmouth  Bass 

317 

Micropterus  salmoides 

Largemouth  Bass 

319 

Pomoxis  niqromaculatus 

Black  Crappie 

331 

Perca  flavescens 

Yellow  Perch 

334 

Stizostedion  vitreum  vitreum 

Walleye  (Yellow  Pickerel) 

338 

Etheostoma  exile 

Iowa  Darter 

339 

Etheostoma  flabellare 

Fantail  Darter 

341 

Etheostoma  nigrum 

Johnny  Darter 

342 

Percina  caprodes 

Logperch 

381 

Cottus  bairdi 

Mottled  Sculpin 
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Table  5:  Sampling  date,  effort  and  total  crayfish  catch  for  the  survey  lakes 


Watershed 

Sampling 
Date 

Species* 

Number  of  Crayfish  Caught 
Females       Males          Total 

Number 
of  Traps 

CUE 

Total 

Lake  Name 

Lake  Code 

CUE 

2CF  Georgian  Bay  Island  Tributaries 

Clearwater                 CL 

Lohi                        LI 

Hannah                    HH 

Middle                     ME 

940824 
940824 
940323 
940823 

Cb 
Cb 

67 
43 
0 
0 

36 
36 
0 
0 

103 
79 

0 

0 

54 
54 
54 
54 

1.91 
1.46 
0.00 
0.00 

1.91 
1.46 
000 

0.00 

2EA  Georgian  Bay  Tributaries 


Grundy 

GY 

920805 

Op 

2 

3 

5 

54 

0.09 

0.22 

Ov 

2 

5 

7 

0.13 

Gut 

GT 

920805 

Cr 
Op 

6 
0 

18 

1 

24 

1 

54 

0.44 
0.02 

0.87 

Ov 

8 

14 

22 

0.41 

Bells  (Silver) 

BL 

920731 

Oo 

18 

15 

33 

53 

0.62 

0.62 

Crawford  (Otter) 

CD 

920730 

Ov 

5 

3 

8 

54 

0.15 

0.46 

Oo 

8 

9 

17 

0.31 

Hams 

HS 

920807 

Cr 
Op 

2 

4 

2 

1 

4 
5 

54 

0.07 
0.09 

0.16 

Old  Mans 

OM 

920731 

Cb 
Op 

1 

1 

2 

1 

3 
2 

54 

0.06 
0.04 

0.66 

Oo 

3 

5 

8 

0.15 

Oi 

7 

15 

22 

0.41 

Simmons 

SS 

920730 

Ov 
Oo 

6 

1 

8 
2 

14 
3 

54 

0.26 
0.06 

0.32 

Black 

BK 

920804 

Cr 

1 

0 

1 

54 

0.02 

0.17 

Op 

4 

2 

6 

0.11 

Oo 

0 

2 

2 

0.04 

Whalley 

WY 

920714 

Op 

1 

0 

1 

54 

0.02 

0.24 

Ov 

7 

5 

12 

0.22 

Bat 

BT 

920807 

Op 

10 

11 

21 

54 

0.39 

0.39 

Rock  Island 

Rl 

920807 

Cr 
Op 

2 

1 

2 
0 

4 

1 

54 

0.07 
0.02 

0.09 

Round  1  (Burton) 

RD 

920804 

Cb 

10 

10 

20 

54 

0.37 

1.31 

Cr 

6 

6 

12 

0.22 

Op 

18 

21 

39 

0.72 

Back 

BA 

920723 

Op 
Ov 

1 
1 

5 
0 

6 
1 

54 

0.11 
0.02 

0.13 

Bell 

BE 

920724 

Op 
Ov 

1 
1 

0 
0 

1 
1 

54 

0.02 
0.02 

0.04 

Fairholme 

FH 

920728 

Cr 

1 

0 

1 

54 

0.02 

0.02 

Long  (Strathdee) 

LG 

920724 

Op 
Ov 

0 

1 

1 
0 

1 
1 

54 

0.02 
0.02 

0.04 

Richmond 

RN 

920723 

Op 

1 

0 

1 

54 

0.02 

002 

Limestone 

LE 

920728 

Op 

2 

1 

3 

52 

0.06 

0.16 

Ov 

4 

0 

4 

0.08 

Oi 

1 

0 

1 

0.02 

*  Cb  -  Cambarus  bartoni: 

Cr  -  Cambarus  robustus; 

Oo  -  Orconectes 

orooinauus: 

Ov  -  Orconectes  virilis 

Oi  -  Orconectes  immunis,    Oo 

Orconectes  obscurus;    Or 

-  Orconectes  rusticus 
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Table  5.  (cont'd) 


Watershed 

Sampling    Species* 
Date 

Number  of  Crayfish  Caught 
Females       Males           Total 

Number 
of  Traps 

CUE 

Total 

Lake  Name 

Lake  Code 

CUE 

2EA  Georgian  Bay  Tributaries 
Fowke  (Spring)  FE 


Manitouwaba 

MA 

Mutton 

MN 

Cheer 

CH 

Frank  ■ 

FR 

Clam  (Big  Clam) 

CM 

Grass  (Sweny) 

GS 

Himbury  (Rock) 

HM 

Kernick  (Pike) 

KK 

920813 

Cr 

40 

Op 

2 

Ov 

1 

Oi 

1 

890816 

Cb 

0 

Cr 

0 

Op 

1 

920721 

Ov 

6 

920813 

— 

0 

910822 

Cb 

1 

Oi 

12 

910828 

Cr 

7 

Op 

8 

910322 

Cb 

13 

910822 

Oo 

6 

910828 

Ov 

3 

80 

6 

1 

1 

1 

1 

3 

8 

0 

2 

26 
17 
13 
28 
16 

5 


54 

148 

1.63 

0.11 

0.02 

0.02 

60 

002 
0.02 
0.05 

009 

53 

0.15 

0.15 

54 

0.00 

0.00 

54 

0.04 
048 

0.52 

54 

0.31 
0.24 

0.55 

54 

0.52 

0.52 

54 

0.30 

0.30 

53 

0.09 

0.09 

2EB  Moon  River  and  Go  Home  River 


Code 

CE 

890815 

— 

0 

0 

0 

20 

0.00 

0.00 

Fogal 

FL 

890810 

Cr 

1 

4 

5 

60 

0.08 

0.55 

Ov 

11 

17 

28 

0.47 

Healey  1  (Conger) 

HY 

890815 

— 

0 

0 

0 

54 

0.00 

0.00 

Oldfield 

OD  - 

890810 

Cr 

2 

9 

11 

60 

0.18 

0.60 

Op 

11 

14 

25 

042 

Armishaw 

AW 

890821 

— 

0 

0 

0 

59 

0.00 

0.00 

Bigwind 

BW 

910813 

Cb 

2 

1 

3 

54 

0.06 

4.27 

Op 

40 

165 

205 

3.80 

Ov 

8 

14 

22 

0.41 

940712 

Op 

209 

251 

460 

53 

868 

994 

Ov 

30 

37 

67 

1.26 

Bird 

BD 

910813 

Op 
Ov 

0 
1 

1 
0 

1 
1 

54 

0.02 
0.02 

0.04 

Chub 

CB 

910709 

Ov 

1 

1 

2 

54 

0.04 

0.04 

940707 

— 

0 

0 

0 

54 

0.00 

0.00 

Dickie 

DE 

910826 

— 

0 

0 

0 

54 

0.00 

0.00 

940708 

— 

0 

0 

0 

54 

0.00 

0.00 

Echo 

EO 

910726 

— 

0 

0 

0 

54 

0.00 

0.00 

Grandview 

GV 

910726 

Cb 

0 

3 

3 

54 

006 

3.13 

Op 

144 

22 

166 

3.07 

IHealey  2  (Macaulay) 

HE 

910820 

Op 

1 

0 

1 

48 

0.02 

0.02 

Leech 

LC 

910820 

— 

0 

0 

0 

54 

0  00 

0.00 

McKay 

MY 

910820 

— 

0 

0 

0 

54 

000 

0.00 

*  Cb  -  Cambarus  bartoni: 

Cr  -  Cambarus  robustus: 

Op- 

Orconectes 

propinauus; 

Ov  -  Orconectes  virilis 

Oi  -  Orconectes  immunis 

;    Oo- 

Orconectes  obscu 

rus; 

Or  -  Orconectes  rusticus 
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Table  5.  (cont'd) 


Watershed 

Sampling    Species" 
Date 

Number  of  Crayfish  Caught 
Females       Males          Total 

Number 
of  Traps 

CUE 

Total 

Lake  Name 

Lake  Code 

CUE 

2EB  Moon  River  and  Go  Home  River 


Moot 

MT 

910809 

Cb 

2 

5 

7 

54 

0.13 

0.13 

940722 

Cb 

2 

4 

6 

54 

0.11 

0.11 

Ril 

RL 

910808 

— 

0 

0 

0 

54 

0.00 

0.00 

Clayton 

CY 

940713 

Cb 

195 

495 

690 

54 

12.78 

12.78 

Axe 

AX 

900726 

— 

0 

0 

0 

54 

0.00 

0.00 

Bittern 

BN 

900802 

Op 

8 

7 

15 

53 

0.28 

0.58 

Ov 

6 

10 

16 

0.30 

Buck 

BC 

900731 

Op 

1 

0 

1 

52 

0.02 

0.37 

Ov 

10 

8 

18 

0.35 

Compass 

CS 

900808 

— 

0 

0 

0 

54 

0.00 

0.00 

Fawn  1 

FN 

900801 

Op 

1 

7 

8 

52 

0.15 

1.13 

Ov 

27 

24 

51 

0.98 

Fawn  2 

FA 

940722 

Ov 

1 

1 

2 

54 

0.04 

0.04 

Fox 

FX 

900809 

Op 
Ov 

0 
2 

2 

1 

2 
3 

54 

0.04 
0.06 

0.10 

Golden  City 

GC 

900726 

— 

0 

0 

0 

54 

0.00 

0.00 

Groves 

GS 

900816 

Cb 

1 

0 

1 

54 

0.02 

0.02 

Haller 

HA 

900803 

— 

0 

0 

0 

54 

0.00 

0.00 

Langford 

LD 

900803 

Op 

1 

1 

2 

53 

0.04 

0.23 

Ov 

5 

5 

10 

0.19 

Lynx 

LX 

900727 

— 

0 

0 

0 

54 

0.00 

0.00 

Rose 

RS 

900816 

Cb 

1 

2 

3 

53 

0.06 

0.06 

Round  2  (McMurrich) 

RO 

900808 

Cb 

0 

1 

1 

54 

0.02 

0.13 

Op 

2 

2 

4 

0.07 

Ov 

2 

0 

2 

0.04 

2EC  Sevem  River 


Bentshoe 

BS 

890719 

Ov 

2 

1 

3 

60 

0.05 

0.05 

Big  East 

Bl 

890717 

Cb 

1 

10 

11 

50 

0.22 

2.80 

Ov 

57 

72 

129 

2.58 

Big  Orillia 

BO 

910730 

Cb 

6 

6 

12 

54 

0.22 

2.59 

Op 

14 

16 

30 

0.56 

Ov 

36 

62 

98 

1.81 

Cinder 

CI 

890809 

Cb 

10 

8 

18 

57 

0.32 

3.43 

Ov 

41 

136 

177 

3.11 

Clear  (Oakley) 

CA 

890810 

Cb 

0 

1 

1 

18 

0.06 

2.89 

Ov 

10 

41 

51 

2.83 

Grindstone 

GE 

910711 

— 

0 

0 

0 

54 

0.00 

0.00 

Gullfeather 

GR 

890721 

Ov 

4 

5 

9 

60 

0.15 

0.15 

Little  Margaret 

LM 

910814 

Cb 

1 

4 

5 

54 

0.09 

3.40 

Ov 

45 

134 

179 

3.31 

Little  Orillia 

LO 

910730 

Op 

4 

2 

6 

54 

0.11 

1.02 

Ov 

26 

23 

49 

0.91 

Nehemiah 

NH 

910814 

— 

0 

0 

0 

54 

000 

0.00 

*  Cb  -  Cambarus  bartoni: 

Cr  -  Cambarus  robustus; 

Op- 

Orconectes 

orooinauus; 

Ov  -  Orconectes  virilis 

Oi  -  Orconectes  immunis 

;    Oo 

-  Orconectes  obscunjs; 

Or  -  Orconectes  rusticus 
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Table  5.  (cont'd) 


Watershed 

Sampling    Species* 
Date 

Number  of  Crayfish  Caught 
Females       Males          Total 

Number 
of  Traps 

CUE 

Total 

Lake  Name 

Lake  Code 

CUE 

2EC  Severn  River 
Red  Chalk  East 


Red  Chalk  Main 


Ridout 
Saw 

Shoe 


RE 


RM 


RT 
SW 

SH 


2HF  Cameron  Lake  Drainage 


89071 1 


940629 


940629 


890818 


Op 

91 

77 

168 

60 

2.80 

6.22 

Ov 

68 

137 

205 

3.42 

Op 

50 

61 

111 

54 

2.06 

4.37 

Ov 

72 

53 

125 

2.31 

Cb 

1 

0 

1 

57 

0.02 

3.81 

Op 

58 

42 

100 

1.75 

Ov 

36 

80 

116 

2.04 

Cb 

15 

9 

24 

54 

0.44 

1.90 

Op 

12 

7 

19 

0.35 

Ov 

37 

23 

60 

1.11 

Cb 

1 

4 

5 

54 

0.09 

0.55 

Op 

9 

16 

25 

0  46 

Cb 

0 

2 

2 

60 

0.03 

0.12 

Op 

1 

3 

4 

0.07 

Ov 

1 

0 

1 

0.02 

Cb 

10 

9 

19 

54 

0.35 

0.35 

Balsam 

BM 

940727 

Cr 

1 

1 

2 

54 

0.04 

0.46 

Op 

0 

3 

3 

0.06 

Ov 

9 

6 

15 

0.28 

Or 

1 

4 

5 

0.09 

Koshlong 

KG 

940805 

Cr 

38 

61 

99 

54 

1.83 

2.15 

Ov 

4 

13 

17 

0.31 

Big  Hawk 

BH 

900718 

— 

0 

0 

0 

54 

0.00 

0.00 

Boshkung 

BG 

890821 

— 

0 

0 

0 

60 

0.00 

0.00 

Goodwin 

GN 

900724 

Cb 
Ov 

1 
2 

0 
1 

1 
3 

54 

0.02 
0.06 

008 

Johnson 

JN 

900720 

— 

0 

0 

0 

54 

0.00 

0.00 

Kabakwa 

KA 

900717 

— 

0 

0 

0 

59 

0.00 

0.00 

Kelly 

KY 

900719 

— 

0 

0 

0 

54 

0.00 

0.00 

Kushog 

KU 

890821 

— 

0 

0 

0 

50 

0.00 

0.00 

Lipsy 

LY 

900815 

Cb 

4 

8 

12 

54 

0.22 

061 

Op 

14 

7 

21 

0.39 

Little  Buckhom 

LB 

900724 

Cb 

21 

38 

59 

54 

1  09 

1.46 

Ov 

8 

12 

20 

0.37 

Little  Hawk 

LH 

900720 

Cb 

1 

1 

2 

54 

0.04 

0.17 

Cr 

2 

2 

4 

007 

Ov 

1 

2 

3 

006 

Little  Kennisis 

LK 

900723 

— 

0 

0 

0 

54 

0.00 

0.00 

Lower  Welch 

LW 

900717 

— 

0 

0 

0 

60 

000 

0.00 

Margaret 

MG 

910717 

Ov 

4 

8 

12 

54 

0.22 

0.22 

Red  Pine 

RP 

900814 

— 

0 

0 

0 

52 

000 

0.00 

Sherborne 

SE 

900717 

Cb 

1 

0 

1 

54 

002 

0.02 

St.  Nora 

SN 

890818 

Cb 

1 

0 

1 

100 

0.01 

0.01 

*  Cb  -  Cambarus  bartoni: 

Cr  -  Cambarus  robustus: 

Op- 

Orconectes 

proDinauus; 

Ov  -  Orconectes  virilis 

Oi  -  Orconectes  immunis:    Oo 

-  Orconectes  obscu 

rus: 

Or  -  Orconectes  rusticus 
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Table  5.  (cont'd) 


Watershed 

Sampling    Species* 
Date 

Number  of  Crayfish  Caught 
Females       Males           Total 

Number 
of  Traps 

CUE 

Total 

Lake  Name 

Lake  Code 

CUE 

2HF  Cameron  Lake  Drainage 
Twelve  Mile  TM 


Basshaunt 


BU 


900718 

— 

0 

0 

0 

72 

0.00 

0.00 

940707 

Op 

0 

1 

1 

54 

0.02 

0.06 

Ov 

0 

2 

2 

0.04 

940805 

Oi 

11 

50 

61 

53 

1.15 

3  36 

Ov 

25 

92 

117 

2.21 

2HG  Scugog  River 
Scugog 


SG 


2KD  Upper  Madawaska  River 


940726 


54 


0.20 


Timberwolf                 TF                940803           Cb              18 

23 

41                54             0.76 

0.76 

Louisa                     LA                940729           Cb              15 

11 

26                54             0.48 

2.55 

Ov             49 

63 

112                               2.07 

*  Cb  -  Cambarus  barton!;    Cr  -  Cambarus  robustus:    Op  -  Orconectes 

proDinquus; 

Ov  -  Orconectes  virilis; 

Oi  -  Orconectes  immunis:    Oo  -  Orconectes  obscurus;    Or  -  Orconectes  rusticus 
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Table  6:  Proportions  of  male  and  female  crayfish  caught  each  year 


Species 

Year* 

Number 

Proportion(%) 

Total 

Species 

Year' 

Number 

Proporti 

on(%) 

_  Total 

of  Lakes 

Female 

Male 

Number 

of  Lakes 

Female 

Male 

Number 

C.  bartoni 

1989 

7 

37 

63 

35 

0. 

oroDlnauus 

1989 

5 

54 

46 

300 

1990 

8 

38 

62 

80 

1990 

7 

51 

49 

53 

1991 

8 

44 

56 

79 

1991 

8 

50 

50 

422 

1992 

2 

69 

31 

16 

1992 

15 

48 

52 

99 

1994 

8 

37 

63 

974 

1994 

7 

45 

55 

619 

Overall 

45 

55 

1184 

Overall 

50 

50 

1493 

C.  robustus 

1989 

3 

18 

82 

17 

0. 

virilis 

1989 

10 

32 

68 

729 

1990 

1 

50 

50 

4 

1990 

9 

54 

46 

123 

1991 

1 

41 

59 

17 

1991 

9 

34 

66 

368 

1992 

7 

46 

54 

126 

1992 

11 

53 

47 

79 

1994 

2 

39 

61 

101 

1994 

12 

39 

61 

774 

Overall 

39 

61 

265 

Overall 

42 

58 

2073 

O.  immunis 

1989 

0 



_ 

0 

0. 

obscurus 

1989 

0 

_ 

_ 

0 

1990 

0 

— 

— 

0 

1990 

0 

— 

— 

0 

1991 

1 

46 

54 

26 

1991 

1 

38 

62 

16 

1992 

3 

38 

62 

24 

1992 

5 

48 

52 

63 

1994 

1 

22 

78 

61 

1994 

0 

— 

— 

0 

Overall 

35 

65 

111 

Overall 

43 

57 

79 

O.  rusticus 

1989 
1990 
1991 
1992 
1994 
Overall 

0 
0 
0 
0 

1 

20 
20 

80 
80 

0 
0 
0 
0 
5 
5 

*  no  lakes  v/ere  surveyed  in  1993 
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Table  7:  Averagesizeof  the  crayfish  species  caught  by  watershed 


species 

Watershed* 

Carapace  Length  (mm) 

Mean 

Range 

Female 

Male  Fl 

Male  Fll 

C,  bartoni 

2CF 

26.4 

26.6 

27.4 

15.0-37.1 

2EA 

32.1 

33.9 

28.0 

16.1  -40.4 

2EB 

34.6 

34.8 

35.3 

18.1-42.7 

2EC 

34.7 

34.6 

30.6 

22.8  -  42.2 

2HF 

28.0 

27.4 

27.6 

13.5-37.7 

C.  robustus 

2EA 

40.0 

42.4 

40.7 

18.1  -54.7 

2EB 

38.7 

42.1 

42.9 

36.2  -  48.2 

2HF 

37.1 

44.9 

40.5 

21.2-50.0 

O.  oroDinquus 

2EA 

24,0 

26.9 

24.1 

17.5-37.0 

2EB 

25.7 

28.6 

244 

11.5-39.5 

2EC 

24.9 

27.1 

23.0 

15.1  -38.0 

2HF 

22.7 

28.1 

25.5 

17.8-34.1 

0.  virilis 

2EA 

30.2 

36.6 

28.8 

17.1  -44.0 

2EB 

31.2 

36.3 

30.9 

21.0-42.0 

2EC 

28.8 

33.3 

28.8 

26.6  -  36.8 

2HF 

27.4 

35.5 

28.5 

21.4-40.6 

2HG 

42.2 

— 

42.4 

14.1  -50.1 

2KD 

24.1 

27.6 

24.7 

17.4-30.9 

0.  immunis 

2EA 

33.4 

33.2 

28.2 

22.2  -  39.8 

2HF 

26.9 

28.7 

26.2 

14.9-33.4 

0.  obscurus 

2HF 

28.3 

35.6 

28.1 

15.2-46.1 

0.  rusticus 

2HF 

35.1 

36.5 

- 

30.9-43.4 

*  see  Table  1  for  watershed  name 
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Table  8:  Crayfish  catch  per  unit  effort  by  species 


Watershed 

Sampling 

Lake  Name 

Lake  Code 

Date 

2CF  Georgian  Bay  Island  Tributaries 

Cleanwater               CL          940824 

Lohi                      LI           940824 

Hannah                 HH         940823 

Middle                   ME          940823 

Crayfish  Species* 


_02_ 


Total 
CUE 


1.91         —  —  —  —  —  —  1.91 

1.46         —  —  —  —  —  —  1.46 

___  —  —  __  0.00 

__  —  —  _  —  —  000 


2EA  Georgian  Bay  Tributaries 

Grundy                   GY  920805  _          _          _          _        0.09         —        0.13  0.22 

Gut                      GT  920805  —        0.43         —          —        0.02         —        0.41  0.86 

Bells  (Silver)               BL  920731  _          _          _        0.63         _          _          _  0.63 

Crawford  (Otter)            CD  920730  _          _          _        0.30         —          —        0.15  0.45 

Hanis                    HS  920807  —        0.07         —          —        0.09         —          —  0.16 

Old  Mans                 OM  920731  0.04         —        0.33       0.13       0.04         —          —  0.54 

Simmons                 SS  920730  _          _          _        0.06         —          —        0.26  0.32 

Black                    BK  920805  —        0.02         —        0.04       0.11         —          —  0.17 

Whalley                  WY  920714  —          _          _          _        0.02         —        0.22  0.24 

Bat                       BT  920807  _           _           —           —        0.39         —           —  0.39 

Rock  Island                Rl  920807  —        0.07         —          —        0.02         —          —  0.09 

Round  1                    RD  920805  0.37       0.22          —           —        0.72          —           —  1.31 

Back                     BA  920723  —          —          —          —        0.11         —        0.02  0.13 

Bell                      BE  920724  _          _          _          _        0.02         —        0.02  0.04 

Fairholme                 FH  920728  0.02         ______  0.02 

Long  (Strathdee)           LG  920724  _          _          _          _        0.02         —        0.02  0.04 

Richmond                 RN  920723  _          _          _          _        0.02         —          —  0.02 

Limestone                 LE  920728  —          —        0.02         —        0.06         —        0.08  0.16 

Fowke  (Spring)             FE  920811  —        1.48       0.02         —        0.11         —        0.02  1.63 

Manitouwaba              MA  890816  0.03         _          —          _        0.05         —          —  0.08 

Mutton                   MN  920721  —          —          —          —          —          —        0.15  0.15 

Cheer                   OH  920811  _______  O.OO 

Frank                    FR  910822  0.04         —        0.48         _          _          —          —  0.52 

Clam                    CM  910828  —        0.31         —          —        0.24         —          —  0.55 

Grass                    GS  910822  0.52         ______  0.52 

Himbury                  HM  910822  _          _          _        0.30         _          _          _  0.30 

Kemick                   KK  910828  ______        0.09  0.09 

2EB  Moon  River  and  Go  Home  River 

Code                    CE  890815  _______  0.00 

Fogal                     PL  890810  —        0.08         _          _          _          _        0.47  0.55 

Healey  1  (Conger)          HY  890815  _______  0.00 

Oldfield                   OD  890810  —        0.18         —          —        0.42         —        0.17  0.77 

Amiishaw                AW  890726  _______  0.00 

Bigwind                  BW  910813  0.05         _          _          _        3.80         —        0.41  4.26 

940712  _           _           _           _        8.68         —        1.26  9.94 

Bird                      BD  910813  —          —          —          —        0.02         —        0.02  0.04 

Chub                    CB  910709  ______        0.04  0.04 

940707  _______  0.00 

Dickie                    DE  910726  ___          —          _          —          _  0.00 

940708  ____  —           _—  0.00 
Echo                     EO  910808  ___          —          _          —          _  0.00 

Grandview                GV  910726  0.06         _          _          _        3.07         —          —  3.13 

Healey  2  (Macaulay)         HE  910820  _          _          _          _        0.02         —          —  0.02 

Leech                    LC  910820  ____          —          _—  0.00 

McKay                   MY  910820  ___          _          ___  0.00 

Moot                     MT  910809  0.13         ______  0.13 

940722  0.11          —           —           —           —           —           —  0.11 

Ril                        RL  910808  ___           —           _           —           —  000 

Clayton                   CY  940713  12.78        ___          _          __  12.78 

Axe                     AX  900726  _______  0.00 

Bittern                   BN  900802  _          _          _          _        0.28         —        0.30  0.58 

Buck                      BC  900731  _           _           _           _        0.02          —        0.35  0.37 

Compass CS  900808  —          —          __          —          _—  0.00 

•  Cb  -  Cambarus  barloni:    Cr  -  Cambarus  robustus:    Op  -  Orconectes  propinguus:    Ov  -  Orconectes  virilis: 
01  -  Orconectes  immunis:    Oo  -  Orconectes  obscunjs:    Or  -  Orconectes  rusticus 
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Table  8.  (cont'd) 


Watershed 

Sampling 
Lake  Code       Date 

Crayfish  Species* 

Total 

Lake  Name 

Cb 

Cr 

Oi          Oo         Op 

Or 

Ov 

CUE 

2EB  Moon  River  and  Go  Home  River 


Fawn  1 

FN 

Fawm2 

FA 

Fox 

FX 

Golden  City 

GC 

Groves 

OS 

Haller 

HA 

Langford 

LD 

Lynx 

LX 

Rose 

RS 

Round  2 

RO 

2EC  Severn  River 

Sentshoe 

BS 

Big  East 

Bl 

Big  Orillia 

BO 

Cinder 

CI 

Clear(Oakley) 

CA 

Grindstone 

GE 

Gullfeather 

GR 

Little  Margaret 

LM 

Little  Orillia 

LO 

Nehemiah 

NH 

Red  Chalk  East 

RE 

Red  Chalk  Main 

RM 

Ridout 

RT 

Saw 

SW 

Shoe 

SH 

900801  _    _    _    _   0.15    —    0.96  1.11 

940722  ______   0.04  0.04 

900809  _    _    _    _   0.04    —   0.06  0.10 

900726  _______  0.00 

900816  0.02    ______  0.02 

900803  _______  0.00 

900803  _    _    _    _   0.04    —   0.19  0.23 

900727  _______  0.00 

900816  0.06    ______  0.06 

900808  0.02    _    _    _   0.07    _   0.04  0.13 

890719  ______   0.05  0.05 

890717  0.22    _____   2.58  2.80 

910730  0.22    _    _    _   0.56    _   1.81  2.59 

890809  0.32    _____   3.11  3.43 

890810  0.06    _____   2.83  2.89 

910711  _______  0.00 

890721  ______   0.15  0.15 

910814  0.09    _____   3.28  3.37 

910730  _    —    —    —   0.11    _   0.91  1.02 

910814  _______  0.00 

890711  _    —    _    _   2.80    _   3.42  6.22 

940629  _    _    _    _   2.06    _   2.31  4.37 

890711  0.02    _    —    —    1.75    _   2.04  3.81 

940629  0.44    _    _    _   0.35    —    1.11  1.90 

940809  0.09    _    _    _   0.46    _   4.35  4.90 

890818  0.03    _    _    _   0.07    _   0.02  0.12 

910704  _______  0.00 


2HF  Cameron  Lake  Drainage 


Balsam 

BM 

Koshlong 

KG 

Big  Hawk 

BH 

Boshkung 

BG 

Goodwin 

GN 

Johnson 

JN 

Kabakvra 

KA 

Kelly 

KY 

Kushog 

KU 

Lipsy 

LY 

Little  Buckhom 

LB 

Little  Hawk 

LH 

Little  Kennlsis 

LK 

Lower  Welch 

LW 

Margaret 

MG 

Red  Pine 

RP 

Shertxjme 

SE 

St.  Nora 

SN 

Twelve  Mile 

TM 

Basshaunt 


BU 


940727  _        0.04         _          _        0.06       0.09       0.28  0.47 

940805  —        1.83         _          _          _          _        0.31  2.14 

900718  _______  0.00 

890825  _______  0.00 

900724  0.02         _____        0.06  0.08 

900724  _______  0.00 

900717  _______  0.00 

900719  _______  0.00 

890821  _______  0.00 

900815  0.22         _           _          _        0.39          _           _  0.61 

900724  1.09         —           _          _          —           _        0.37  1.46 

900720  0.04       0.07         _  _          _           _        o.04  0.15 
900723  _______  0.00 

900717  _______  0.00 

910717  ______        0.26  0.26 

900814  _______  0.00 

900718  0.02         ______  0.02 

890818  0.01         ______  0.01 

900718  _______  0.00 

940707  _          _          _          _        0.02         _        0.04  0.06 

940805  _           _        1.15         _           _           _        2.21  3.36 


2HG  Scugog  River 

Scugog  SG         940726         ______        0.20       0.20 

2K0  Upper  Madawaska  River 

Timbenwolf  TF  940803       0.76         ______        0.76 

Lousia LA  940729       0.48         _____        2.07       2.55 

*  Cb  -  Cambarus  bartoni:    Cr  -  Cambarus  robustus:    Op  -  Orconectes  propinquus:    Ov  -  Orconectes  virilis: 
Oi  -  Orconectes  immunis:    Oo  -  Orconectes  oljscurus;    Or  -  Orconectes  rusticus 
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Table  14:  One-way  ANOVAs  separating  lakes  based  on 
C.  bartoni  presence-absence 


C.  bartoni  Presence/Absence 

One-way  Analysis 

of  Variance 

Variable 

Model  Sum    Model  Error  Mean 

Error 

Model 

F 

P 

of  Squares       DF         Square 

DF 

r-squared 

Value 

Value 

Elevation 

3009 

0.308 

98 

0.091 

9.760 

0.002 

Lake  Area 

0.097 

0.098 

98 

0.010 

0.991 

0.322 

Lake  Volume 

0.048 

0.236 

98 

0.002 

0.203 

0.653 

Maximum  Depth 

0.902 

0.135 

98 

0.064 

6.681 

0011 

Alkalinity 

0.635 

0.260 

98 

0.024 

2445 

0.121 

Al 

0.002 

0.003 

98 

0.007 

0.719 

0.399 

Ca 

0.358 

0.173 

98 

0.021 

2.069 

0.154 

CI 

0.028 

0.015 

98 

0.019 

1.919 

0.169 

Cond25 

0.104 

0.032 

98 

0.032 

3.288 

0.073 

Die 

0.117 

0.042 

98 

0.028 

2.773 

0.099 

DOC 

0.001 

0.000 

98 

0.041 

4.150 

0.044 

Fe 

0.157 

0.034 

98 

0.044 

4.554 

0.035 

K 

0.013 

0.008 

98 

0.015 

1.535 

0.218 

Mg 

0.347 

0.155 

98 

0.022 

2.235 

0.138 

Mn 

0.770 

0.597 

98 

0.013 

1.289 

0.259 

Na 

0.088 

0.040 

98 

0.022 

2.192 

0.142 

NH4 

0.083 

0.725 

98 

0.001 

0.115 

0.736 

N03 

0.374 

1.839 

98 

0.002 

0.203 

0.653 

TKN 

0.002 

0.000 

98 

0.035 

3.601 

0.061 

pH 

0121 

0.039 

98 

0.031 

3.131 

0.080 

S04 

0.089 

6.871 

98 

0.000 

0.013 

0.910 

TP 

0.099 

0.026 

98 

0.037 

3.799 

0.054 
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Figure  5 
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Figure  6 
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